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Abstract: Primordial non-Gaussianity is generated by interactions of the inflaton field, 
either self- interactions or couplings to other sectors. These two physically different mech- 
anisms can lead to nearly indistinguishable bispectra of the equilateral type, but generate 
distinct patterns in the relative scaling of higher order moments. We illustrate these classes 
in a simple effective field theory framework where the flatness of the inflaton potential is 
protected by a softly broken shift symmetry. Since the distinctive difference between the 
two classes of interactions is the scaling of the moments, we investigate the implications for 
observables that depend on the series of moments. We obtain analytic expressions for the 
Minkowski functionals and the halo mass function for an arbitrary structure of moments, 
and use these to demonstrate how different classes of interactions might be distinguished 
observationally. Our analysis casts light on a number of theoretical issues, in particular we 
clarify the difference between the physics that keeps the distribution of fluctuations nearly 
Gaussian, and the physics that keeps the calculation under control. 
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1. Introduction 

The next few years will bring an extraordinary amount of new information about the pri- 
mordial cosmological fluctuations. The Planck satellite is measuring the temperature 
anisotropies of the Cosmic Microwave Background (CMB) with substantial improvements 
in precision, reaching down to smaller length scales. Current and upcoming large volume 
surveys are cataloging cosmological structures (eg - galaxies, clusters and voids) whose 
number and distribution depend on the same initial conditions that source the CMB tem- 
perature anisotropies |2|, |3|, |4], [5L [BL [7L pL Pj. These primordial fluctuations provide a re- 
markable observational window into the history of the early universe and offer the potential 
to probe fundamental physics at the energy scale of inflation, orders of magnitude above 
what can be reached in the laboratory. To make the most of this data, it is important to 
understand in detail what microscopic information is encoded in the primordial curvature 
fluctuations, and how this may be extracted from an observational perspective. 

The increasing observational precision will allow us to significantly constrain, or even 
measure, statistics of the primordial fluctuations beyond the power spectrum, collectively 
called non-Gaussianity. To date, theoretical work and observational constraints of non- 
Gaussianity have focused mostly on signatures of the 3-point correlation function (the 
bispectrum), which provides a leading indicator of non-Gaussian effects in most scenarios. 
Higher order moments are challenging to compute and even more challenging to constrain 
observationally. But, when the primordial curvature fluctuations are calculated beginning 
with a particle physics model for the inflaton field(s), and especially when non-Gaussianity 
is observably large, there are relationships between the moments of the fluctuations coming 
from the structure of the original theory. This is an especially striking feature of the most 
robust and appealing inflationary models: they rely on some symmetry that, although 
weakly broken, determines much of the structure of the terms appearing in the theory. 
When non-Gaussianity is large, this structure can be directly mapped to the pattern of 
primordial fluctuations in the gravitational potential and the order-by-order perturbative 
expansion contains evidence of the underlying symmetry. More than the presence or ab- 
sence of a three-point function, it will be consistency relations between the moments them- 
selves (and ideally with the evolution of the inflationary background) that will make a more 
compelling case for inflation or force us to consider alternatives. 

In this paper we discuss two types of consistency relations that are typical of particle 
physics inflationary scenarios. One is a feature of nearly every non-Gaussian model consid- 
ered in the literature so far and the other is quite new. We work, for the most part, within 
one of the simplest and most popular microscopic frameworks: inflation driven by a single 
scalar field with an underlying shift symmetry. Although the starting point is simple, there 
is a rich array of possible non-Gaussian signatures. Focusing on global properties of the 
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correlation functions rather than just the bispectrum, we uncover two distinctive classes of 
interactions: 

• When non-Gaussianity comes from self-interactions of a single field, the correla- 
tion functions are typically hierarchical: the n-th dimensionless moment scales like 

/ l/2\ n-2 

A4 n oc ( IV ^ J where is the amplitude of the fluctuations and X measures the 
strength of the interaction (X ~ Jnl)- 

• When non-Gaussianity is generated by interactions with other sectors, the dimen- 
sionless moments exhibit a very different structure. They may fall off much more 
slowly, or not at all. This novel scaling of the moments indicates that such scenarios 
are more non- Gaussian than their hierarchical counterparts, given an amplitude of 
the three-point function. 

We demonstrate how the two classes may be easily confused at the level of the three-point 
function, but would be particularly distinguished by observables sensitive to the difference 
in the pattern of moments, for example the statistics of collapsed objects. 

We also clarify the difference between the physics that keeps the Probability Distri- 
bution Function (PDF) nearly Gaussian, and the physics that keeps the calculation under 
control. In the hierarchical case, the structure of self-interactions is intimately linked to 
the validity of perturbation theory for the interactions (ie, keeping loop corrections to the 
two-point small). In contrast, the case with non-hierarchical scaling can be under control 
even in a regime where the dimensionless moments are not only non-hierarchical but are 
not even a converging series. This case is important because it helps us to understand the 
theory of inflation, although the simplest models are already ruled out observationally. 

Although we will calculate features of a particular class of models, our goal is to 
illuminate where the programme of constraining or detecting non-Gaussianity is headed 
more generally. For example, a detection of non-Gaussianity from a particular bispectral 
template would be fantastic, but would not uniquely determine the primordial physics. We 
would then be faced with the need to do at least one of several difficult measurements: a 
very precise measurement of the bispectral shape, a measurement of the four point function, 
or a measurement that would fold in information about higher order moments. The first 
two are technically challenging but straightforward, so it is the last option that we will 
explore in this paper. 

We will find that our analysis is also relevant for a more general issue: any constraint 
or detection of non-Gaussianity should be confirmed with observables that are sensitive to 
non-Gaussianity but that are not direct n-point function measurements. For example, we 
might cross-check measurements of Minkowski functionals in the CMB or the distribution 
of massive objects in the late universe (the mass function). To make a useful comparison 
between different statistics, we must extract the parameter measuring the strength of the 
interaction. Absent a complete theoretical prediction for all the moments, this means in 
practice that we expand about a Gaussian in the higher moments. Any valid truncation 
of this expansion depends on the structure of the moments, so we must be cautious about 
interpreting the results when a particular expansion has been assumed. 
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As a case in point, it is tempting to try to use rare objects, the most massive clusters, 
to constrain or find evidence for non-Gaussianity [10, 11, [l^, 13, 14|. Such measurements 



provide valuable information on scales different from CMB observations and, moreover, 
they are subject to a different set of systematics. While there are many statistical issues 
in making such an analysis convincing [15, [If], [l?], [18| , we will show here that we must also 



be very careful in interpreting what such objects imply for the physics of the primordial 
fluctuations. On the other hand, given good enough observations, the structure of the 
cumulants could be extracted from Large Scale Structure (LSS) measurements, yielding a 
remarkable window into the interaction structure of the underlying field theory. 

This paper is organized as follows. In section ^ we review single field inflation with 
an approximate shift symmetry, using effective field theory to classify the possible non- 
Gaussian effects. Section |3] discusses the scaling of the correlation functions that is associ- 
ated with each type of non-Gaussianity. Section Q considers another example of what we 
will call the feeder mechanism. Next, Section ^ considers the observational implications 
and discusses the possibility to disentangle various classes of microscopic interactions. We 
present generalized expressions for the halo mass function and for the Minkowski func- 
tionals, suitable for scenarios with non-hierarchical scaling. Finally, in section ||, we con- 
clude. Appendix A reviews the classification of inflationary models, Appendix B provides 
a detailed computation of the non-Gaussian correlation functions from inverse decay and 
Appendix C provides some details on the computation of non-Gaussian correlators from 
rescattering. 

2. Inflation with a Shift Symmetry 

The standard "vanilla" inflationary models contain only a single scalar field with La- 
grangian 

C = -l( dip f- V (cp) + --- (2.1) 

where • • • denotes other fields, including the Standard Model of particle physics. In order 
to support a long quasi de Sitter phase with \H\ <C H 2 (here H = a/a is the Hubble scale 
and a(t) is the scale factor) we require a flat scalar potential, quantified by the slow roll 
parameters 

Ml {V'\ 2 V" 

< = ^{v) ■ i = M *T- (2 - 2) 

In addition to the classical condensate, (p(t) = (<p(t, x)), quantum fluctuations are inevitably 
present due to the inflationary background: 

<p(t,x) = <t>(t) + 5ip(t,x.) . (2.3) 

The (super-horizon) curvature perturbation is obtained by a gauge transformation 

C(t,x) = -?Mi,x). (2.4) 
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Thus, in this scenario the same field supports the background evolution and also generates 
the observed primordial curvature fluctuations. The power spectrum, V^k), is defined by 



(C 2 (x)) = J ^P C (k). (2.5) 

1/2 H 2 

In the vanilla scenario we have TV (k) = * where quantities on the right-hand-side 
are understood to be evaluated at horizon crossing, k = aH. Prom measurements of 
the amplitude of fluctuations in the cosmic microwave background | |19| |, ~ 2.42 x 



10 U002MPC- 1 J with the s P ectral index n s ~ 1 ~ °- 034 - 

In this simple framework we can identify two important scales: H and \J~eM v (the 



scale appearing in the conversion (|2.4j ) between Sip and £). These scales characterize most 
of the interesting low energy physics associated with the scalar curvature fluctuations. In 
addition, the Planck scale, M p & 2.4 • 10 18 GeV, controls quantum gravity effects and 
provides a natural UV cut-off. 

In the simple model (|2.l| ), the only interactions are gravitational or else arise through 
the slow roll potential. In both cases, the interaction strength is controlled by slow roll 
parameters and non-Gaussianity is too small to be observable in the next decade or more 
|2C, ^l], 22]. Additional interactions, controlled by a new scale, / <C M p , will arise once we 



try to embed the theory (2.1) into a more complete framework. 



2.1 Effective Field Theory of the Inflaton 



The slow roll parameters (2.2) are notoriously sensitive to Ultra- Violet (UV) physics and 



e, \r]\ <C 1 may require significant fine tuning. An appealing way to achieve slow roll in a 
naturally is to invoke a softly broken shift symmetry 

ip — > p + const . (2-6) 



In the simplest models <f)(t) can be considered as a clock and ( |2.q ) leads to the approximate 
time-translation invariance that one might regard, in the context of the full gravitational 
system, as the "true" symmetry of inflation. 

From a field theoretical perspective, the most cogent way to implement the symmetry 
(|2.G| ) is to assume that the inflaton is a Pseudo-Nambu-Goldstone-Boson (PNGB) which 
results from the spontaneous breaking of some global symmetry in the matter sector, at 
the scale /. The low energy theory will be characterized by the Goldstone mode, ip, 
along with any other light fields which are present, for example gauge fields. The original 
implementation of this idea was Natural Inflation p3| , 24, ^5|, 26], and many variations 
on this theme have followed @ |7|, [§, ||| [H], |32], ||, @, ||, g |7|, g|, g|, [I| ||, g2|, 
motivated in part by the fact that such models may naturally generate an observably large 
amplitude of primordial gravitational waves. 1 Despite this significant body of work, it was 



observable primordial gravitational wave background requires that the inflaton moves a super- 
Planckian distance in field space during inflation plf . Ensuring that the potential is flat over such large 
displacements seems unnatural, unless there is extra physics at work such as an underlying shift symmetry. 
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only recently that the cosmological perturbations in slow roll models with an underlying 



shift symmetry were understood comprehensively [44|. 

Following the standard logic of effective field theory, we write down terms in the action 
which are consistent with the symmetry (|2.6|): 



Co = ~\{drf - ±<pG»& - ±<pFF + f>^£^ + (2.7) 

Here F^ v = d^Ay — dyA^ is the field strength for some U(l) gauge field A^, while G a ^ v is 
the field strength associated with some non-Abelian gauge group. (The gauge fields are 
assumed to have standard kinetic terms which we omit for ease of presentation.) 

The Langrangian ( |2.7| ) should be understood as an effective theory below the scale /. 
At higher energies, the global symmetry is restored and one must include additional degrees 
of freedom beyond the PNGB. Generically, the symmetry breaking scale / also controls 
the strength of interactions in (2.7); absent fine tuning (or the appearance of additional 



scales in the problem) we expect c n = 0(1) and a = 0(1). 

Notice that the shift symmetry (^^) strongly constrains the allowed interactions, and 
so will have important implications for observable non-Gaussian signatures. Goldstone's 
theorem implies that <p can only have derivative interactions. 2 These include the self- 
interaction terms (dip) n and also the pseudo-scalar couplings to gauge fields, ipFF and 
tpG a G a , that we have already included in fl2.7|) . 3 In addition to these leading order terms, 
higher derivative couplings to gauge fields are also allowed, suppressed by more powers of 
/. Moreover, couplings to fermions will also be generically present; these will play no role 
in our discussion. 

To drive inflation, we must slightly break the shift symmetry (|2.6[), This may be done 



in a controlled way, so that the smallness of the slow roll parameters (2.2) is protected by 



the mildness of the symmetry breaking effects. There are two approaches: 

1. Non-Perturbative Breaking: Instanton configurations of the non-Abelian gauge 
field will generically break the shift symmetry, even if it is exact at the classical level. 
The coupling ipG a G a leads to a potential of the form 

V np (<p) = A 4 cosf^ +••• (2.8) 



where A is a nonperturbatively generated scale and • • • denotes higher harmonics. 
Notice that a discrete subgroup of the original symmetry remains unbroken: cp — > 



Explicit breaking: One may instead break the symmetry ( |2.6D explicitly at tree 
level. We parametrize such effects by a non-periodic potential V ex ((p) which may, 



for example, have a power-law form yii 4 p tpP. See [06| , 37 1 for an interesting field 



theory construction. Explicit breaking may also be motived from a string theoretic 



perspective, for example, by wrapping branes on suitable cycles [32 



See j46|] for a review. 

3 The quantity FF is a total derivative so that ipFF gives a derivative interaction d^K^ after integrating 
by parts. Similarly, also, for the non-Abelian gauge field. 
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Putting everything together, including symmetry breaking terms, we are led to consider 
the following Lagrangian: 

C = ~\{drf - V(ip) - JjL<pFF + £ c„^L_ + . . . ( 2 .9) 

^ n=l •* 



F(^) = y ex (V9)+A 4 cos(^)+--- (2.10) 



where the potential is 

V(in\ = V„Jin\ 4- A 4 rns ( ^ 

Note that in this simple framework there is a single scale, /, controlling the interactions 
in ( |2.9[) . The introduction of explicit symmetry breaking effects (via V ex ) may introduce 
an additional scale. However, explicit symmetry breaking effects are assumed to be para- 



metrically small, in order to preserve the slow roll conditions (2.2). Hence, we expect that 
any additional interactions induced by explicit symmetry breaking will be too weak to 
contribute an observable non-Gaussian signature. 

2.2 Computational Control 

In the last subsection, we have argued that Q2,9|) provides a generic low energy effective 



description of a PNGB inflaton, neglecting gravity. This theory is characterized by a 
single scale, /. However once we couple ( |2.9[ ) to gravity, new scales arise in the problem 
associated with the curvature (Hubble) scale, its rate of evolution, and the Planck mass. 
In the inflationary context we should re-examine what it means to call the complete theory 
"effective". 

We have in mind that ip is the Goldstone mode associated with a global symmetry 
broken at the scale /. In this case, we should presumably require / <C M p since otherwise 
symmetry breaking would occur above the quantum gravity scale where Quantum Field 
Theory (QFT) is expected to become invalid. Moreover, it is expected that quantum 
gravity effects (such as the virtual appearance of black holes) will explicitly break the 
original global symmetry [^] and lead to corrections to fl2.9Q that are suppressed by powers 
of f /M p . Note that the condition / <C M p seems also to be a requirement of string theory 
p7[| , and is consistent with the "gravity as the weakest force" conjecture pSf . 

As discussed previously, we can identify two important scales - H and ^feM p - which 
characterize most of the physics associated with the curvature perturbations in the vanilla 
scenario. The non-gravitational sector ( p.Sj ) contains an additional scale, /, which arises 
naturally once we try to embed single-field inflation into a more complete UV framework, 
so it is helpful to organize our discussion using the ratio of the scale associated with field 
theory interactions, /, to the scales associated with gravity, H and y/eM p . We introduce 
the dimensionless quantities 

PM P = — J — , PH = -J ■ (2.11) 

These are of course not both independent of the amplitude of fluctuations, but they are 
useful to organize the physics. 



4 See P5|] for a proof of this statement. 
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In order for the effective description (|2.9| ) to make sense, we should require that inflation 
takes place below the scale of the symmetry breaking that gives rise to the PNGB inflaton. 
This implies that fin *C 1. For larger values of H we would need to include additional 
degrees of freedom, for example the massive field whose Vacuum Expectation Value (VEV) 
breaks the original global symmetry. 

Control over the low-energy description (^^) often also requires that the higher di- 
mension operators fall off with increasing powers of some small parameter. Otherwise, we 
would need to re-sum an infinite series of such terms, the explicit form of which is not, in 
general, known. The condition (3h <C 1, which we have already discussed, helps to control 
higher derivative interactions involving the fluctuations, {d8ip) n . This is easily seen by esti- 
mating the derivatives as § ~ y ~ A more stringent constraint on model parameters 
comes from requiring that we can truncate the derivative expansion in (^1]) when studying 
the dynamics of the homogeneous background, <f)(t). In this case one has 



The higher derivative terms in (2.9) modify the sound speed, c s , of the inflaton fluctuations. 



Note that the condition ( 2.12 ) implies that the sound speed is not too different from unity: 



1-4 (^) 2 (f) 2 = 1-4&p&«1. (2.13) 



It is possible that some underlying symmetry principle strongly constrains the structure 
of derivative self-couplings, making it possible to re-sum the series and obtain c s < 1 in a 



controllable way. DBI inflation [49] provides an example of such a theory, where the UV 
symmetry is associated with higher-dimensional Lorentz invariance |50t| . One might also 
invoke Galileon symmetry |pl| . 

In summary, we generically expect that the low-energy effective description fl2.9|) will 
be under computational control provided (3h *C 1 and PhPm p <C 1. Notice that there is a 
priori no restriction on the size of (3m p by itself, however, we do require / <C M p . In most 
interesting cases, the conditions for the validity of the effective field theory ( |2.9| ) can be 
summarized as: 

H < / < M p . (2.14) 
2.3 Mechanisms for Non-Gaussianity 

Although our low energy effective description ( |2.9[ ) is simple, it nevertheless allows for a rich 
array of non-Gaussian signatures. Note that the same shift symmetry (|2.6|) which protects 



the flatness of the potential ( 2.10 ), also strongly constrains the interactions of the inflaton. 
Hence, this symmetry has important implications for non-Gaussianity. There are three 
kinds of interactions which may be distinguished, each leading to a unique non-Gaussian 
effect: 



5 The latter condition can be relaxed in cases where we have more information about the UV structure 
of the theory. 
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Derivative Self-Interactions: The higher derivative self-interactions (dip) n lead 
to an equilateral bispectrum with nonlinearity parameter Jnl^ c~ 2 p^ j. An obser- 
vationally interesting signal is possible for c s <C 1. 

Resonant Effects: A different source of non-Gaussianity arises due to the periodic 
contribution to the potential (|2.10| ), which has its origin in nonperturbative symme- 
try breaking effects. The self-interactions encoded in this term give rise to a novel 



oscillatory signal in the bispectrum 53, 54, [)5|, 56] 



Pseudo-Scalar Couplings (Inverse Decay, Feeder field): The pseudo-scalar 
interaction with Abelian gauge fields, (pFF, leads to a nonperturbative production 
of gauge quanta which, in turn, source non-Gaussian inflaton fluctuations via inverse 



decay [44, 57]. This mechanism was first understood in [44] and is reviewed in Ap- 
pendix B. Non-Gaussianity is observable when / <C M p , which seems required in any 
case for the consistency of the effective field theory description. This scenario is a 
realization what we will call the "feeder" mechanism, discussed in detail below. 

In all cases the interactions which source non-Gaussianity occur within (or close to) 
the horizon. Moreover, in all cases there is no significant iso-curvature perturbation on 
large scales, nor super-horizon evolution of £. 

2.4 A Note About "Single Field" Inflation 

Before moving on, we pause to define how we distinguish between single field and multi- 
field models. Obviously, it is not useful to define a model as "single field" when only one 
degree of freedom appears in the Lagrangian: any model that incorporates realistic particle 
physics must necessarily include many spectator fields which play no role in the inflationary 
dynamics. A more cogent definition is to classify a model as "single field" whenever there 
is only one field driving the background inflationary dynamics and that same field is also 
responsible for generating the super-horizon curvature perturbation, £. Throughout this 
paper it is understood that the expression "single field" has this meaning. Appendix A 
elaborates on this definition, establishing a taxonomy of inflationary models that we find 
useful to organize the relevant physics for the background and fluctuations. This is also 
relevant for understanding where the "feeder" mechanism sits in relation to the "single 



field" effective field theory of the fluctuations [58]. 

Notice that there is no obstruction to allowing a "single field" inflaton to couple directly 
to matter fields, as would seem to be a pre-requisite for successful reheating |59|, |60|, |6l], |62| . 
(Here we use the word "matter" to denote any particles that are not directly involved in 
driving the background inflationary dynamics.) Hence, the model (|2.9|) is classified as 



single field, despite the fact that we have explicitly included couplings to gauge fields. 
2.5 The Feeder Mechanism: Unexpected Results in Single Field Inflation 



The pseudo-scalar coupling discussed in subsection 2.3 is actually an example of a larger 



class of scenarios that are "single field" according to our definition, but which exhibit novel 
and unexpected phenomenology [44, 57, |63|, 64, |65|, |66]j . In particular, these models may 
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give rise to observably large non-Gaussianity without recourse to any of the phenomena 
that add a layer of explanation on top of the inflationary shift-symmetry, including effects 
such as small sound speed [j52fl , higher derivatives |37|, |6£j] , special initial conditions [52 



|69[ |, potentials with sharp features |70], [52]], dissipative effects [71], fine-tuned inflationary 
trajectories, post-inflationary effects (such as preheating [[7^, |73|, [74|), etc. 

The key mechanism that underlies this novel class of single-field inflation models works 
as follows. Prom an effective field theory perspective, we generically expect direct couplings 
between the inflaton and one or more matter fields, as would seem to be required for 
successful reheating. In some cases, the coupling to the time-dependent condensate can 
excite quanta of the matter field. These produced matter particles may, in turn, source 
inflaton fluctuations via the same coupling. Inflaton fluctuations that are generated in this 
way tend to be highly non-Gaussian, and are complementary to the usual (nearly) Gaussian 
quantum fluctuations from the vacuum. When this effect is observationally significant, we 
refer to the matter field as a feeder, since its fluctuations grow in the environment of the 
quasi de Sitter background and, in turn, feed the inflaton perturbations, helping them to 
grow large, non-Gaussian and strong. 

To date, several models have been explored which exhibit the feeder mechanism ||44| , 
57, [35], |66], 75]. 6 These models differ in the details of the particle production (either 

tachyonic instability E[ or an instantaneous violation of adiabaticity |6^]) and also the 
feed-back mechanism (inverse decay [44] or rescattering Q). However, the basic physics 



is very much analogous. We expect that many other examples can be found. In [75, 76] 
the production of gravitational waves via the feeder mechanism was studied. 

In all known examples, the "feeding" dynamics take place entirely on scales inside or 
comparable to the horizon and the bispectrum is closer to equilateral than local. There is no 
significant iso-curvature perturbation on super-horizon scales. For k <C aH the curvature 
perturbation depends only on the inflaton field and is frozen, as in the vanilla scenario, so 
the feeder models really are "single field" in the conventional sense, both at the level of 
the homogeneous background dynamics and also at the level of the perturbations. 

3. Structure of Correlation Functions 



In subsection |2.3| we discussed three mechanisms - sound speed, resonance and inverse 
decay - which can give rise to observable non-Gaussian effects in the context of single field 
inflation with an underlying shift symmetry. In this section, we will study the structure of 
correlation functions from each of these mechanisms, showing how this pattern of moments 
encodes information about the underlying microphysics. In section ||, we will discuss how 
this information might be extracted from an observational perspective. 

3.1 The Dimensionless Moments 

To compare models we will work with the amplitude of the dimensionless moments, strictly 
defined from the equilateral limit of scale-invariant momentum space ra-point functions. We 



5 See Appendix B for a review of the model proposed in 
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write the n-th cumulant 7 (the connected part of the real space n-point function) as 

d 3 ki d 3 k 2 d 3 k n 



(C"(X)) C EE 



(2vr) 3 (2tt) 3 (2tt) 3 



(2 7 r) 3 P n , c (A; J ) ( 5( 3 ) ^ki 



vi=l 



(3.1) 



When the underlying fluctuations are exactly scale invariant, the polyspectra scale as 
Pni ~ A;" 3 ^" 1 ). The dimensionless moments, M n , are defined by 

_ k^-VP^jk, fc, ■ ■ ■ , fc) _ k^-VPn^k, k,...,k) 

[^(fc)]"/ 2 " (27T^f/ 2 ' ^ 

where is defined by (|2.5| ) and we have taken the equilateral limit | k\ | = | | = • • • = 
|fc n | = fc. For a Gaussian distribution M. n >z = 0, so these provide a measure of non- 
Gaussianity. 

Note that because the M n are dimensionless, we can use them as the dimensionless mo- 
ments for either the primordial curvature £ or for the matter era Bardeen potential, $ = |£ 
by replacing P n £ — > P n ^ and — > A|. Throughout, we will ignore any scale-dependence 
(including n s / 1) that will make Ai n a function of scale. Such scale dependence is typi- 
cally mild, although it may still be relevant for understanding the perturbation theory if 
it is stronger at higher n (eg, see fT^| ) and is likely to be observable if the non-Gaussianity 



is not too small 78, [7^, |80| . Note that, for resonant models, this convention amounts 
to ignoring the oscillatory functions that modulate the amplitude. In Section || we will 
move the discussion to real space, and then these approximations can be avoided, but a 
smoothing scale is typically introduced which also changes the scaling slightly. Since the 
Ai n defined above establish the scaling of the coefficient of some momentum integrals even 
in the real space case, we will often schematically write A4 n as without specifying 

the argument of £. 

In this paper we are interested in how the pattern of dimensionless moments encodes 
the physics of the underlying inflaton interactions. A useful diagnostic of the pattern of 
moments is the ratio 

^±1 (3.3) 

We will focus on the parametric dependence of this ratio on the small parameter that 
controls the strength of interactions, showing that very different patterns arise for self- 
interactions as compared to interactions with other sectors. The ratio (3.3) will also depend 



on numerical coefficients which may grow with n due to combinatorics. Such model- 
dependent characteristics are, of course, different from the parametric pattern of moments 
and must be considered on a case-by-case basis. In practice, observations of rare objects 
are only sensitive to relatively low order moments (certainly n < 10 is sufficient), so we do 
not expect combinatorics to modify our key results in any significant way (see section |H|). 
We leave a more detailed analysis to future works. 

We will call a set of moments {A^ n } ordered if M n+ \ < Ai n for all n. Often, higher 
order moments fall off with increasing powers of some small parameter and so are naturally 



7 Here we have omitted any smoothing, which is irrelevant since our discussion will focus only on the 
scaling of the integrand, for the time being. In section Hour definition of the moments will be made rigorous. 
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ordered. In this case we can develop a systematic expansion of the Probability Distribution 
Function (PDF) around the Gaussian; see section |5[ Moreover, non-Gaussian corrections 
to the PDF are controlled by the smallness of the dimensionless moments. 



3.2 Self-Interactions and Hierarchical Scaling 



Our effective field theory description (|2.9| ) contains two distinct classes of interactions: self- 
interactions (leading to resonance or sound speed effects) and couplings to other sectors 
(leading to inverse decay effects). In the simple scenario we consider, both kinds of effects 
are controlled by a single scale, /. Shortly, we will verify explicitly that the the self- 
interaction terms in ( |2.9| ) lead to a specific pattern of moments referred to as hierarchical. 
To shed light on the underlying physics, it is interesting to first consider this point from a 
somewhat more general perspective. 

Suppose the fundamental Lagrangian contains only interactions that couple £ to itself 
and, moreover, there is a single scale defining the interaction strength, I in appropriate 
units. Heuristically, we expect the ratio of the (n + l)-function to the n-function to scale 
as 

(C n+1 ) 

, A . /c oc 1 ■ 2tt 2 V c . (3.4) 

(C n )c C V ' 

The factor of I is obvious: higher correlation functions involve either higher order inter- 
action terms or additional vertices, which necessarily carry more powers of the interaction 
strength. The factor of arises because the (n + l)-th correlator has an additional ex- 
ternal line, as compared to the n-th correlator. To connect this external line to the rest 
of the diagram requires an additional propagator, which scales as Q 2 ~ Vr- We define our 
normalization so that the 2-point function gives one power of Vt and the 3-point defines 
one power of the interaction strength. Putting these together, we expect 

(O oc (If- 2 (2TT 2 V ( ) n ~ 1 . (3.5) 



Then the dimensionless moments (3.2) scale as 



n — 2 



M n oc (X 2 2tt 2 P c ) 2 . (3.6) 



This is the scaling commonly referred to in the literature as hierarchical, and the moments 
are ordered provided 

X • V2irVl /2 < 1 . (3.7) 

In many examples, this same condition also controls the loop expansion. 

If there are additional scales in the problem, or fine-tuning, it may be possible to obtain 
non-hierarchical moments from self-interactions. Nevertheless, as we will see explicitly, the 
heuristic argument which we have presented does apply to the most widely-studied non- 
Gaussian mechanisms. 
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3.3 The Local Ansatz 

In order to illustrate the statements made in the last subsection, let us consider the so-called 
local ansatz for the curvature perturbation: 

( = C 9 + \fNL [C 2 g ~ <0] , (3.8) 

where Q g is a Gaussian random field. This expression provides a useful phenomenological 
description of non-Gaussianity and may also arise in multi-field or curvaton models. The 
cumulants scale as (C n )c ~ {Inl)™ (C^) n 1 an d hence the dimensionless moments have the 
structure 

n-2 
2 

(3.9) 

where A n is some numerical coefficient from combinatorics. The local ansatz reproduces 
( |3.6| ) with I ~ /nl- The moments are ordered when ( |3.7D is satisfied. Equivalently, we 
can write fNLiCg) 1 ^ 2 ^ 1- This shows that the naive condition for the non-Gaussian part 
of ([T^) to give a subdominant contribution to the two-point function also ensures that the 
moments are ordered. 

In general, scenarios with hierarchical scaling need not give a non-Gaussianity of the 
type (|3.§|), Rather, the bispectrum may have a shape which is very different from local. 
Very different bispectra are often compared by defining an effective nonlinearity parameter 
as 

f eS ^ (C ) /o 1n \ 

~ 3/£2\2 ' ( 3 - 10 ) 

where the bispectrum has been evaluated in the equilateral limit and the factor of 5/3 makes 
this convention consistent with the usual local ansatz. For some models, this definition is 
consistent with I ~ fff L so that our generic result, equation (|3.6[ ), is brought into the 
familiar form 

(3.11) 

However, we will see in two examples below (the inflaton/curvaton scenario and reso- 
nant non-Gaussianity) the definition of fff L directly from the bispectrum does not always 
correctly capture the interaction strength, while the ratio of consecutive dimensionless 
moments (|3.3| ) does. It is worth emphasizing that the role of fff L in our discussion is 
to provide a simple measure of the parameter which controls the strength of interactions. 
The full bispectrum contains information beyond the amplitude - the shape and running of 
non-Gaussianity are both important discriminating tools - but these are not the subject of 
our attention here. Indeed, we will ultimately be interested in observables such as the halo 
mass function and Minkowski functionals, which are relatively insensitive to the detailed 
shape of non-Gaussianity. 

3.4 Inflaton/Curvaton non-Gaussianity 

An interesting variation on the local model comes when two fields contribute to the cur- 
vature, one of which has local-type non-Gaussianity. In other words, an Antatz of the 



M 



local 



3 f 

TJNL 

5 



2tt 2 P c 



OC 



\JNL, 
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type: 



C(x) = (j)(x) + a(x) + \f NL o(x) 2 (3.12) 
5 



Defining the ratio of power in the two-point from the two fields as 



.2 V a V a 



f = - ° = -f . (3.13) 



we see that the leading terms in the higher moments go like 

n-2 

These moments have the hierarchical scaling, although with an extra parametric depen- 
dence on £ 2 in each of the dimensionless moments so that defining fff L as in Eq.( 3.1C| ) 



above would be misleading as a measure of the interaction strength. Notice that if the 
real space expansion Eq.( |3.12| ) had no term linear in a above this would become the "un- 



Gaussiton" which scales like the feeder mechanism we discuss beginning in subsection 3.8 
below. 

3.5 Sound Speed Effects 

Let us now consider a less trivial example: non-Gaussianity generated by the derivative 
self-interactions in fl2.9|) . Models of this type have been well-studied in the literature and 
it is known that the strength of interactions is controlled by the sound speed, c s , and that 
large non-Gaussianity is possible when c s < 1. (In terms of the scale / suppressing the 
derivative interactions, small c s indicates / much closer to H than to \/eM p [ ]81~| , |82|.) The 
dimensionless moments for generic derivative self-interactions were estimated in [[77], |8T| 

71 — 2 

M% = A n (c7 4 2tt 2 ^) 2 , (3.15) 



where A n is independent of model parameters. This is exactly the form (3.6) with I 



Fnl ~ c s~ 2 ' consistent with the expectation that the sound speed should be associated with 
the strength of interactions. 

If we consider an action with an arbitrary function of powers of the field and its first 
derivative, we would find additional parameters that describe the amplitude of other non- 
Gaussianity terms at each order n. However, each parameter generates a family of higher 
moments with the hierarchical scaling. Explicit examples can be found in p3| ]. If the 
structure of higher derivative terms is altered by imposing an additional symmetry, the 



scaling may also change [33|. 
3.6 Resonance Effects 



The periodic contribution to the potential ( 2.10| ), which arises due to nonperturbative sym- 



metry breaking effects, provides another interesting source of non-Gaussianity. Expanding 
the field in terms of fluctuations as (I2.3T) we have interaction terms of the form 



C n = -p^{8vY (3.16) 
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where the oscillatory time dependence of the effective coupling is encoded by 

(-l) n / 2 cos (W) n even; 

w " 3 < -.1)^1.7^ -odd. (3 ' I7) 



Non-Gaussianity arises due to a resonance between the oscillation frequency of the effective 
coupling and the wavenumber A: of a given mode. Resonance occurs inside the horizon, 
leading to a bispectrum that is closer to equilateral than local. 

The higher order correlation functions due to the interaction ( |3.16| ) were recently com- 
puted |5j| . At leading order the result is 

MT- = A n /3~J p (/4 p 2vr 2 P c ) V , (3.18) 

where A n are dimensionless numbers and {3m p was defined in ( ^.ll[) . 8 Hence, we recover 
the hierarchical structure ( j3lj| ) with X = 0\t = (\feM p / f) 2 . 

It is interesting, also, to compare Pi) with (|3~15|) . From Eq.fljTl|), it looks like 
Pjj plays the role of 1/c 2 in the small sound speed models. However, there is an extra 
dependence on the parameter (3m p in the coefficient of each M. n , so that the amplitude of 
the three point is given parametrically by instead of fi\ Ip - That is, 

<C 3 >«<(2vr 2 P c ) 2 . (3.19) 

As with the inflaton/curvaton example, this highlights the fact that amplitude of the third 
moment less to distinguish models than information about the structure of the moments 
and their parametric dependence on the scales in the problem. 9 To check or uncover the 



pattern in Eq.( 3.18 ) one would need at least measurements of the three and four point 
amplitudes. 

Note that our discussion of resonance effects has focused on Fourier space amplitudes. 
Here we characterize the size of non-Gaussian moments by the amplitude of their oscillatory 
polyspectra. Care may be required to translate this particular discussion into real space 
where integration of the oscillatory functions over k can lead to cancellations. 

A final interesting point about the resonant non-Gaussianity is that at high enough 
order n, the correlation functions may be dominated by a contribution that scales slightly 
differently, with one less power of (3m p - When (3m p is a large number, this improves the 



convergence of the series of moments [55]. 
3.7 Inverse Decay Effects 

So far, we have seen that the self-interactions in (|2.9| ) give rise to a hierarchical structure of 
correlation functions where the dimensionless moments fall off according to Eq.([3,6|). This 



8 The quantity /3m p = y/eMp/ f was called a in reference |5q ], 

9 The full bispectrum contains more information than the skewness; the shape and running of non- 
Gaussianity are also important discriminators. However, for the kinds of probes that we will ultimately be 
interested in - the abundance of rare objects and morphological indicators - it is the integrated moments 
that are most relevant. More on this later. 
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exhausts two of the three mechanisms which may generate observable non-Gaussianity in 
the context of single field inflation with an underlying shift symmetry (resonance and sound 
speed effects). We now explore the structure of correlation functions associated with the 
third mechanism: inverse decay effects. We will see that a dramatically different kind of 
scaling behaviour arises in this case. 

Inverse decay effects generate non-Gaussian correlations in a way that is fundamentally 
different from the kind of self-interaction vertices that were discussed in subsection |3.2| . 
Hence, we do not expect that equation (|3.6| ) will hold. The key physics of inverse decay 
was first understood in [Q and is reviewed in Appendix B where we also compute the 
parametric dependence of all ra-point functions. Let us briefly re-capitulate the key aspects 
in order to provide a heuristic estimate of the cumulants. 

We are interested in the Lagrangian ( |2.9| ), assuming that resonance and sound speed 
effects are negligible. The calculation proceeds in two steps. First, consider the coupling 
of the gauge field to the time-dependent condensate 

£D-> 2 -^FF + ... (3.20) 

For (j) = const the second term is a total derivative which has no impact on the classical 
equation of motion. However, with 0^0 this term breaks the conformal invariance of 
the gauge field action at a scale M* ~ ^j-. This leads to a tachyonic production of gauge 

fluctuations with 5A^ ~ He w ^. Here £ = ^hj ~ TT quantifies the scale associated with 
the breaking of conformal invariance as compared to the Hubble scale. 

The next step in the calculation is the feed-back of the produced gauge field fluctuations 
into the inflaton perturbation. The relevant terms in the Lagrangian are 

£ D —{pSipf - - ^StpFF + ■■■ (3.21) 

The last term mediates inverse decay processes where two of the produced gauge fluctua- 
tions generate an inflaton mode 5A^ + 5A^ — > 6<p. Clearly, this term has a very different 
structure from the self-interaction vertex considered in subsection |3.2| . The amplitude of 
this effect can be naively estimated from the equation of motion 



V 2 

d 2 t + 3H 8 t -^+m 2 



5 V = ^FF. (3.22) 



The relevant dynamics takes place near horizon crossing, so we estimate the derivatives as 
d ~ H. This, along with our previous estimate SA^ ~ ife 7 ^, suggests (J^inv.dec ~ ^j—e 2n ^. 

Using C ~ ~jS(f we are led to expect (C"> oc (^^e 27r? )". 

In Appendix B we explicitly compute the parametric dependence of all non-Gaussian 
n-point functions generated by inverse decay. The result for the polyspectra is 

P n ,(k) = (0.0031-^^V (3.23) 



"2ttH / £ 4 J A;3(--i) 
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See equation ( |B-28 ). Here A n are dimensionless numbers, independent of model param- 



eters. By construction A% = \ and it may be verified that A3 = 0(1) [44|. We expect 
that the higher order coefficients are also order unity, although combinatorial factors may 
modify this expectation at very large n. This is an interesting issue, but since most observ- 
ables are not sensitive to extremely high order moments we will leave this matter to future 
investigation. Equation ( [3.23 ) confirms the heuristic estimate in the last paragraph (up to 



a power-law dependence on £ which arises from carefully solving for the time-dependence 
of the gauge field fluctuations). 

In addition to the inflaton fluctuations generated by inverse decay, we also have the 
usual quantum vacuum fluctuations from inflation, <5c^ vac ~ H (the homogeneous part of the 
solution of ( |3,22D ) . These are nearly Gaussian and hence irrelevant for the n = 3 and higher 
moments. However, the 2-point function gets important contributions from both inverse 
decay and also vacuum fluctuations; see equation ( B-27] ) for the explicit result. There are 



two limiting cases which must be distinguished, depending on which effect dominates 

{ ,„ H . n ; „ , vacuum dominated 

2 (3 24) 

(2 ^2^2 2^ (o.0031^^p^ , inverse decay dominated 



where the spectrum is defined as in (|2.5| ). The condition to be in the vacuum dominated 
regime is 

aH e 2 ^ r , , 

0.0031— — — j— <C v27r [vacuum dominated] (3.25) 



which is evident from comparing the two limiting behaviours in ( |3.24| ); see also ( B-27| ). 



(Obviously, the inverse decay dominated regime corresponds to ( p.25| ) with the inequality 
reversed.) 

3.8 A Different Kind of Scaling 

We are finally in a position to compute the structure of dimensionless moments from inverse 
decay effects. It is straightforward to show that 

^inv.dcc = f K (^ir^fr) . vacuum dominated . 
[ A n , inverse decay dominated 

Let us consider each of these limiting behaviours separately. 

1. Vacuum Dominated Regime: Here the condition ( |3.25| ) is satisfied and the mo- 
ments are easily seen to be ordered in the sense that M%Z? ec < AC v - dec . However, 
the structure of correlation functions is very different from what would be obtained 
for self-interactions. To see this, we introduce an effective nonlinearity parameter as 
in equation ( |3.10| ). The dimensionless moments can now be written as 

AC v ' dcc oc l n (3.27) 
oc [f%iV(- j [vacuum dominated] 
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where in the first line we identify X = °^ 31 ®jLs?li as the natural measure of the 



strength of interations in this scenario. In the second line of (|3.28[) we re-write the 
moments in a way that is appropriate for contrasting with small sound speed models, 
Eq. (gHD , and uses Eq.( |3.10| ) to define the effective non-linearity parameter. For a 
given value of fff L the moments from inverse decay fall off much more slowly than for 
self-interactions, so that even if both mechanisms are tuned to give a 3-point function 
of the same size, the inverse decay scenario will lead to a PDF that is intrinsically 
more non-Gaussian. 

2. Inverse Decay Dominated Regime: In this case the dimensionless moments are 
not even ordered] This can be understood physically: when inverse decay effects 
dominate, £ is bi- linear in the Gaussian field SA^. Hence, we do not expect the PDF 
to be close to Gaussian in any meaningful sense. Remarkably, the computation which 
leads to ( |3.26| ) is still under control. In Appendix B we discuss the constraints coming 
from calculability. The most important consistency condition arises from ensuring 
that backreaction effects are negligible. This amounts to the constraint 

r- H 2 aH e 2 ^ 

V2nl = 0.0031—^^-^ « 0(1 • (3-28) 

27T\(f>\ J 4 

This constraint ensures that the n-point functions of £ are small. However, in the 
inverse decay dominated regime the small parameter cancels out of A4 n , leading to 
a non-ordered structure. 

We see, then, that inverse decay effects are remarkably different from the self inter- 
actions discussed in subsection |3.2| . In the latter case, the same small parameter which 
controls the hierarchical structure of moments is also responsible for ensuring computa- 
tional control. In the inverse decay case, on the other hand, these two statements are 
effectively decoupled. It is possible to have non-ordered moments in a computationally re- 
liable regime. We expect that similar results will hold also in other examples of the feeder 
mechanism; more on this in the next section. 

3.9 Salvaging the Inverse Decay Dominated Regime 

The inverse decay dominated regime is theoretically interesting, since it provides a novel 
example of a non-hierarchical PDF which is computationally reliable. However, this sce- 
nario gives /^jf 1 XL 10 3 , which is strongly disfavoured from an observational perspective. 
There are several ways to circumvent this difficulty. One idea, put forward in |85| , is to 
consider M 3> 1 U(l) gauge fields, each coupled to the inflaton as in our model fl2.9| ) with 
the same coupling strength a/ f. Because these gauge fields are statistically independent, 
the n-point correlation functions are multiplied by M. To hold fixed the normalization of 
the power spectrum we must require that the quantity X oc scales as A/"" 1 / 2 . 

Given this relation, the effective nonlinearity parameter must scale as A/"" 1 / 2 . For £ — > oo 
the nonlinearity parameter in the A" = 1 model saturates to fff L — > 8400. Hence, in this 
regime we can write the result for the multi-field case as 

ff L = 8400A'- 1 /2 j (3.29) 
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so that J\f ~ 10 3 is sufficient to render the model compatible with observational bounds on 
CMB scales. The dimensionless cumulants M n scale as AA 1-n / 2 , giving 

M% = (l.2 • lO-Vfi)"" 2 M*= l . (3.30) 

Remembering that Ai^ =1 are independent of model parameters, we see that the cumulants 
exhibit a scaling very similar to the standard hierarchical scaling, Eq. ( |3 .11 ). The key 



difference is that the small parameter V^ 2 fffi L has been replaced with 1.2 • 10 4 /^- 
4. Another Example: Production of Massive Particles 

In section || we studied the structure of correlation functions for the various scenarios that 
can lead to observable non-Gaussianity in the model ( |2.9| ). We found that inverse decay 
effects lead to a radically different structure than self-interactions. We expect that this 
novel pattern of moments will be shared by other examples of the feeder mechanism. To 
verify this, let us now turn our attention to a model that was studied in detail in 63, 64. 65 1 



(see (6(| for a review). This model is somewhat tangential to our original construction of 
a PBNG inflaton ( |2.9| ). Nevertheless, we will see that many similarities arise. 
Consider the theory 

C = -\(drf ~ V&) ~ \(dx? ~ ^{<p - 6>)V , (4.1) 

where (p is the inflaton and x is a matter fields. This theory provides a simple and explicit 
toy model to investigate the implications of particle production during inflation and shares 
many qualitative similarities with the inverse decay mechanism discussed previously. Mod- 
els of the type (|4.1|) are quite natural in the context of open string inflation and various 



extensions of this simple framework have been studied in W& 86, p 



Let us briefly discuss the dynamics of the model (4.1). At the moment when <fr(t) = (fro 
the effective mass of x vanishes, leading to an instantaneous violation of adiabaticity and 
associated particle production. The variance of produced fluctuations is 

g\4> - M 

where the number density of produced quanta is n x ~ a~ 3 (g<p) s ^ 2 where the factor a~ 3 
reflects the usual volume dilution of non-relativitic particles. The dominant effect of par- 
ticle production of the observable spectrum of curvature fluctuations arises because the 
produced, massive x particles can rescatter off the condensate to generate bremsstrahlung 



radiation of long- wavelength 6<p fluctuations 63]. This process is described by the equation 



V 2 

d\ + 3Hd t - -y + m 



5^ = -g 2 [^t)-Mx 2 - (4.3) 



(In Appendix C we discuss the theory of this equation in more detail.) Rescattering leads 
to a localized, bump-like feature in the primordial power spectrum. This may be described 
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by the following simple fitting function [34|: 



VAk) 



Kq 



+ 



\ 3 J 



(4.4) 



where ko is the pivot and kf, is the location of the feature. The first term is associated 
with the usual vacuum fluctuations from inflation (the homogeneous solution of ( |4.3| )) 
while the second term arises due to rescattering processes (the particular solution of (|4.3| ) ) . 
The bump-like feature in the power spectrum ( |4.4[ ) will be associated, also, with localized 
non-Gaussian features in the n-point correlation functions. 



The polyspectra for the model (4.1) are very far from scale invariant, reflecting the 
fact that rescattering effects are most important for modes leaving the horizon near the 
moment when <f> = (f)Q. To compute the cumulants in this model, we smooth on scales 



Rb ~ k b , corresponding to the location of the bump in (4.4). The calculation is carried 



out in detail in Appendix C. Here we simply quote the key result: 



(C(x)) R6 ~ a n 



(2tt 5 ) 



H 2 



M<t>\ 



(4.5) 



where a n are dimensionless coefficients; more on these shortly. The variance of fluctuations 
gets important contributions from both rescattering and also the usual vacuum fluctuations. 
We can distinguish two regimes, depending on which effect dominates: 



<C 2 (x))i4 ~ 



H 



(2ir) 2 <t> 2 

a 2 (2irg) 



2 H 4 



(2tt) 2 4> 2 



vacuum dominated 
rescattering dominated 

The dimensionless moments are defined in position space as Ai n = . 
following result 



(4.6) 



We find the 



Ml 



a n (27rg) 7 



a n (a2 



-n/2 



vacuum dominated 
rescattering dominated 



(4.7) 



We see a structure very much analogous to (|]2|), which was obtained for inverse decay. 
However, unlike the previous case there is an extra dependence on model parameters in 

the prefactor: a n = {jj^j A n , where A n is a numerical coefficient, independent of model 
parameters. Taking this into account, we have an extra parametric dependence in the 

vacuum dominated regime: Ai n oc A n C {2irg) n where C = .In the rescattering 



n-2 
2 



dominated regime, we instead have M. n oc (C -1 ) 2 . In the case at hand C > 1 and so the 
moments are ordered. The situation is reminiscent of ( p.30 ) where the moments fall off at 
a rate similar to the hierarchical case, however, the small parameter is decoupled from the 
amplitude of the variance. 



The regime where rescattering effects dominate the variance (4J3) is under computa- 
tional control. The bump-like feature in ( |4.4D can dominate over the usual scale-invariant 
contribution, even while the coupling is perturbatively small (g 2 < 1) and backreaction 
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effects are negligible Note that nothing forbids the inclusion of a bare mass term fj?x 2 
in the Lagrangian ( |4.1| ); even if the x held is classically massless at <fi = <po one expects 
that such a term will be generated by radiative corrections. The new parameter fi has the 
effect of reducing the efficiency of particle production effects. For models obtained from 
super-gravity or string theory it is natural to have ~ H. In this case the polyspectra are 
essentially unsuppressed and our analysis is not modified in any important way |65| . 

The analysis presented in this section suggests that our qualitative results for the 
pattern of moments may be rather generic for any realization of the feeder mechanism. It 
would be interesting to confirm this intuition explicitly. 



5. Observational Implications of the Structure of Correlation Functions 



The different scaling of the moments above means that the physical origin of non-Gaussianity 
may only be evident by looking at observables beyond the bispectrum. The four-point func- 
tion in the CMB has been studied for nearly a decade [||, |89| and recently constraints at 
an interesting level have been obtained. For local type non-Gaussianity four different tech- 
i 93 1 have been used to constrain the amplitude to a level where the 



9C, 91 



mques 

scaling of the moments might soon be checked. In addition, Fergusson et al [92] apply their 
technique to also constrain other trispectra, including the equilateral type. If we know 
which shapes of trispectra and bispectra to compare from particular theoretical models, 
direct comparison of the amplitudes of those moments is a promising direction, especially 
given the expected sensitivity of the Planck satellite |j89[ . In this section we investigate 
a different approach, demonstrating how cross-checks of Gaussianity like Minkowski func- 
tionals and cluster number counts (which use integrated moments rather than momentum 
space n-point functions) are sensitive to the structure of the correlation functions. 



5.1 Summary of Scalings 

Before considering the smoothed, real-space moments we summarize our findings based on 
the analysis of momentum space amplitudes of the polyspectra. We find two particularly 
important patterns in the dimensionless moments in terms of the parameters I labeling the 
interaction strength and the amplitude of fluctuations V^. One scaling is a characteristic 
of single-field interactions, called the hierarchical scaling. The other is generated by what 
we call the feeder mechanism and is characteristic of interactions between the inflaton and 
a second field that is a spectator in terms of the inflationary background. The scalings are 
captured by the ratio of consecutive dimensionless moments, but the moments themselves 
may have amplitudes shifted by numerical coefficients or parametric coefficients that do 
not scale with order n. These are: 



• Hierarchical scaling (local model, derivative self-interactions): 
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• Hierarchical with parametric coefficient (resonance self-interactions): 



M n oc A n C (1 2 2tt 2 V c ) 2 



• Feeder scaling (inverse decay): 

M n oc A n l 



Feeder with parametric coefficient (massive particle production): 



M n oc A n C V 



We illustrate the first three scalings in Figure |T|, where for visualization purposes param- 
eters have been chosen so that the numerical value of M3 is equal for all scenarios (but 
unrealistically large) and all the A n = 1 for simplicity. 
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Figure 1: A comparison of the scaling behaviors. All cases have been normalized to give the same 
value of M3. The purple diamonds show the standard hierarchical scaling, generated by derivative 
self-interactions and the local ansatz, for example. For a fixed amplitude of the bispectrum, scenarios 
of this type are the least non-Gaussian of the scalings discussed here. The brown circles show a 
hierarchical scaling with an extra parametric dependence in the individual moments, as occurs 
in resonant non-Gaussianity or the mixed inflaton/curvaton scenario. This last case should also 
mimic some of the difference that might come from having very different A n factors. The blue 
stars show the feeder field scaling given for example from the inverse decay, particle production and 
"un-Gaussiton" scenarios. 
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Above we have focused on the simplest scenario where the contribution of the feeder 
field to the power spectrum is sub-dominant. In the opposite regime the moments are 
not generically ordered and a systematic expansion of the PDF around a Gaussian is not 
possible. Moreover, the non-ordered regime is ruled out observationally in the simplest 
models. 

5.2 Mass Function 

Perhaps the most interesting implication of the different scaling of the moments is in 
what we would predict for the non-Gaussian mass function of gravitationally collapsed 
objects. Following the Press-Schechter approach |)4|] (which has previously been extended 
to non-Gaussian initial conditions in |)|, |97|, [7^, ^]), to predict the number density 
dn(M)dM of halos with masses in the range (M, M + dM) we need an expression for the 
non-Gaussian probability distribution of (normalized) density fluctuations smoothed on 

/ „ \l/3 

a scale R associated to halos of mass M by R = I -r—M I . Labeling the distribution 



\ 47r/9 

P(v, M) with v = 5 / an - here a R is the smoothed variance defined, explicitly below - the 
fraction of volume in collapsed objects (halos) is 

roc 

F(M) = / dvP(u,M) (5.1) 
with 5 C the threshold for collapse. The mass function is 

— (M,z) = -2i-^, (5.2) 
dM y M dM y ' 

where p = Sl m Pcrit is the average (co-moving) matter density. The smoothed density field 

is given by 

d 3 k 
(2tt) e 

where Wn(k) is the Fourier transform of a window function, commonly taken to be a 
top-hat in real space giving 

Jxr 3sin(A;i2) 3cos(kR) 

WR{k) = (5-4) 

The relation between the primordial curvature perturbation £ and the linear perturbation 
to the matter density 5 = 5p/p today is 

5(k,z) = M(k,z)((k) 



/ 7^W R (k)5(k,z) (5.3) 



M(k,z) = D(z)T(k)e (5.5) 

where D(z) is the linear growth function, z is the redshift, and T(k) is the transfer function. 
The smoothed variance and 3-point are, for example, 

/d 3 k f d 3 k' -, 
(2^)3 J J^W R (k)W R (k')M(k,z)M(k',z)(C(k)C(k')) 

T W R (kfM(k,z) 2 V ( (k) (5.6) 



o 



{S3r) = Iw?!w?I |^^^^M 1 M 2 M 3 (C(A; i )Cfe)C(A : 3)) 
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To find an approximate expression for the non-Gaussian probability density, we will use 
the Edgeworth expansion. Although one can think of various ways to express a slightly non- 
Gaussian distribution as a series (the Gram-Charlier and Gauss-Hermite expansions are two 
others) the Edgeworth expansion is neatly organized by the moments of the distribution 
and is a true asymptotic expansion, allowing one to calculate the error when the series is 
truncated. The expansion was originally developed to describe distributions that approach 
the Gaussian (according to the central limit theorem) as the number of degrees of freedom, 
n, increases. In this case the expansion is organized in powers of 1/y/n. Although only the 
first few terms are typically quoted in the physics literature, Petrov |9^, |100| (references 
are the English translations) developed an expression for the complete series and proved 
the expansion is asymptotic for cases where a ~ 1 j \fn. 

The common hierarchical scaling often prompts a definition of the "reduced cumulants" 



Sn,R = }l R n-l ( 5 ' 7 ) 
\ 6 R 



2 \n—l 



which are related to the dimensionless moments by 

M^ R = S n , R a n R ~ 2 . (5.8) 

Note that the dimensionless moments are always redshift independent. Then Petrov's 
expression, written in terms of (a generic) a rather than ^/n and in terms of the smoothed, 
reduced cumulants, is 

P Md ^^{ 1+ |^ ( g H , + , M n^(^) fe "} . m 

Here H n (y) are Hermite polynomials defined by H n (v) = (— l) n e^ 2//2 g^-e -iy2//2 and r = 
k\ + A>2 + • • • + k n where the set {k m } is built of all non- negative integer solutions of 

k 1 + 2k 2 -\ \-nk n = s. (5.10) 

We can re- write this expression in a form that will allow us to organize the expansion 
appropriately for non-hierarchical scalings by using the dimensionless smoothed moments 
M n> R. Making use of Eq.( |5.10D , we find 



I s=l {k m } m=l Ky ' I 



— "-"^(l + px{v, R) + P2{v, R) + ... 



2vr 

where pi(is,R) organize the expansion in terms of the small parameter which controls the 
strength of interactions. As we motivated for several examples in earlier sections, it is the 
smallness of the Ai n that controls the departure from Gaussianity, and so to keep good 
control of the errors in the expansion we should group terms of the same order according 
to the scaling of the moments. 
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For the hierarchical scaling, where it is natural to group the M\r term with the 
C^3 b) 2 term, etc, this leads to the standard grouping in powers of a (~ 1/y/n for n 
degrees of freedom) |fT8| . Then the first few terms look like 



p[ h \u,R) = -^H 3 (u) (5.12) 

{h) ( u \ tt I \ , 1 

P 2 \ v i R ) = -^r H ^ u > + 2 




On the other hand, for the feeder field type scaling we expect the {-M-^ R ) 2 term to be 
roughly of the same order as the A4q r term, etc. In that case, we have 

p[ f \u,R) = ^H 3 (u) (5.13) 
p ( 2 f \u,R) = ^j&Htly) 
p { /\v,R) = ^H 5 (u) 

{f), m (m{,rI ( M {r\ \ rrr \ 

p£ (v,R) = I -gf- + 2 [— ) I H *(y) 

m {Ml R Mt R M f 7R \ 

p [ i\^R) = ( tw. ' +^rj H ^) 

We should mention that the integrations smoothing over the momentum space n— 1 spectra 
(as in Eq.(p\^)) appear to change the scaling of the moments slightly (eg, see Figure 1 in 
§H where M± :R » 2.5M 2 3R ). O ne can probably account for this systematically, and in 
any case that contribution to the scaling will be the same for all scenarios, regardless of 
how the coefficients of the momentum space spectra scale. Since both the hierarchical and 
feeder field scenarios would shift in the same way, this should not qualitatively affect our 
comparison of the two cases. 



We note that expansions following from Eq.(5.11) need only resemble that of the actual 
non-Gaussian distribution near the peak, with a smaller range of validity onto the tail 
the more the actual distribution deviates from a Gaussian (see [ |101| ] for some illustrative 
examples). However, the point where the errors are large can be computed by looking at 
the size of the terms in the series. From this we also see that the inverse-decay dominated 
regime, where all the M n are the same size, is not well-captured by any truncation of this 
expansion except very near the peak (y <C 1). 



Using these probability distributions in Eq.(5.2), we can expand the collapse fraction 



F(M) = F {M) + F X (M) + F 2 {M) + ... (5.14) 
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where 



Fo( M) = i E rfc(^) 



(5.15) 



and the non-Gaussian terms F\(M) etc contain dimensionless moments grouped at each 
order according to their scaling. To write the mass function we will need derivatives of 
those terms with respect to mass (or smoothing scale). The first term is the usual Gaussian, 
so regardless of non-Gaussian type we have 



dFh 
dM 



da v c 



2tt dM a 



(5.16) 



Then the ratio of the non-Gaussian Edgeworth mass function to the Gaussian has the same 
structural form for either scaling: 



ting , 



F[ hJ) '{M) F^ J) '(M) 



I Edgeworth 



1 + 



(hj)'/ 



F^M) 



+ 



F>(M) 



+ 



(5.17) 



In the next section we use a specific example to show how the difference in mass 
function might help us distinguish the two scenarios observationally. Although we have 
not explicitly written all the terms above, we will show the results up to third order for for 
the feeder field case to show how the series is behaving. 

For reference, the first few Hermite polynomials are: 



(5.18) 



flo(f) 


= 1 






H x (y) 


= V 






H 2 (v) 


= v 2 


- 1 




H 3 (is) 


= v 3 


- 3^ 




H A (v) 


= u 4 


- Qu 2 + 3 




H 5 (i>) 


= u 5 


- 10z> 3 + 15^ 






= v % 


- lhu A + 45z^ 2 - 


- 15 


H 7 (v) 


= u 7 


- 21i/ 5 + 105i/ 3 


- 105^ 


H 8 (is) 


= u 8 


- 28i/ 6 + 210i/ 4 


- 420^ 2 + 105 


H 9 (v) 


= u 9 


- 36z/ + 378i/ 5 


- 1260z/ + 945^. 



5.3 Comparing Hierarchical Scaling to the Feeder Mechanism 

As a reference point, let's assume that a bispectrum of the equilateral shape has been 
detected. That is, we have measured the amplitude, B eq , using the equilateral template 
|l02j l (for simplicity we assume n s = 1) 



ah)ak 2 )ak 3 : 



i 



(kik 2 



- k 2 + k 3 ) B eq 
+ 2 perm. — 



(5.19) 



(hk 2 k 3 



+ 



k\ 



+ 5 perm. 
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Suppose we find B eq = 3.3 x 10 , corresponding to an effective nonlinearity parameter 
fx L = B eq /(6V^) ~ 100. If we attribute this to a model with hierarchical scaling, then we 
expect M\ ~ 0(M 3 ) 2 . On the other hand, if we attribute it to a feeder field (eg - inverse 
decay) then we expect M.{ ~ 0{M 3 ) 4l I z (ignoring for the moment the slight discrepancy 
in the exact shapes from the two scenarios). Clearly, to distinguish these scenarios we need 
more information. This could come from a measurement of the trispectrum, 10 but here we 
will explore a consistency check from the distribution of high mass/high redshift clusters. 

Since we have measured the three-point in our scenario the first non-Gaussian term 
in the mass function is determined. To distinguish the scenarios we will compare what we 
would predict for the non-Gaussian mass functions at next order. That is, we will compare 
the second order mass function with hierarchical scaling to the second order mass function 
with feeder field type scaling. 

We need the non-Gaussian terms for each type of scaling. For the hierarchical case 



F[ h) {M) 



F?\m) 



e 2 



Ah) 
l R,3 



2^ 3! 



-^-H 3 (u c ) + - 



M 



(h) 

R,3 



3! 



(5.20) 



and for the feeder field 



F[ f) {M) 



F 2 (/) (M) = 



M 



if) 



2tt 3! 
-A *>(/) 

"2w 4! 



'-H 2 (v c ) 



H 3 {v c ) 



(5.21) 



For a scale invariant equilateral shape, the dimensionless moments M.^,Ri M.a,r are only 
very weakly dependent on the smoothing scale. Then the derivatives of the non-Gaussian 
terms simplify: 



5 r dFj da 



a 2 dv dM 

and so the hierarchical and feeder field terms are 



(5.22) 
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Ah)' 
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Af)> 
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R3 e 2 da v r . -vl 



3! J2^ dM a 



H 3 (u c ) = Fq 



M%\ 1 
-t H 'M + 2 



771/ ~^R,3 



M 



(/) 



R,A 



M 



(h) 

R,3 



3! 



(h) 

R,3 



3! 



H 3 (v c 



(5.23) 



H 6 (v c 



3 The precise shape of the bispectrum and the tensor-to-scalar ratio may also be a valuable discriminators. 
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Putting these expressions into Eq.( |5.17| ) gives the Edgeworth prediction for the non- 
Gaussian correction factor to the mass function. Figure |2| shows the ratio of the non- 
Gaussian mass functions to the Gaussian for the two scenarios with the approximation 
where the terms proportional to derivatives of the A4 n ,R are dropped. 

z=0 z=0.5 

dn/dM dn/dM 




Figure 2: The ratio of the non-Gaussian to Gaussian mass function for redshifts z = and z = 0.5. 
Each panel shows a set of lower curves that assume equilateral non-Gaussianity corresponding to 
Inl = 100 and upper curves corresponding to Jml = 250. In each set of curves, the black solid lines 
shows the result for the pdf truncated at M^.r while the blue dashed and purple dotted lines show 
the second order result for the hierarchical and feeder field scalings respectively. Notice that the 
range of masses where the fourth moment is a significant correction is model dependent. Finally, 
the red dot-dashed lines shows the next order {M$,r) correction for the feeder field scenario. 



From Fig. ||, there are two relevant differences between the scenarios. First, in the 
feeder field case the correction from moments above the third one is relevant at significantly 
lower mass. Second, the hierarchical scaling predicts fewer high mass/high redshift objects 
than the non-hierarchical scaling, although at a level that might require futuristic data sets 
to distinguish. 

A template for an equilateral type trispectrum has been constrained from the CMB 
by Fergusson et al [p2[ . The template, chosen to match (up to a numerical factor) one of 
several explicit single-field shapes calculated by Chen et al [ |103| ], is 

(C(ki)C(k2)C(k 3 )C(h)} = (2ir) 3 5 3 D (h + k 2 + h + h) P4,c(*i, fc 3, h) (5.24) 

4 



<C(£i)C(k)C(&)C(feO) = (^fs'Uh + h + h + iU) 8t^(2^ 2 p c ) 3 -L TT I 

A i=X Ki 

K = (k 1 + k 2 + k 3 + /c 4 )/4 



Ref. 1 92] quotes the constraint: 

t NL = (-3.11 ± 7.5) x 10 6 . (5.25) 

From this constraint one can work outs how a particular small sound speed model like DBI 
is constrained compared to the inverse decay scenario. 
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5.4 Log-Edgeworth Expansion 

Although the Edgeworth expansion of the probability distribution has several useful fea- 
tures, it leads to a mass function that is often not well-behaved. Recently LoVerde and 
Smith [98] have proposed a slightly different expansion technique which leads to what 



they call the 'Log-Edgeworth' mass function 11 . They compare this mass function and 
the standard Edgeworth expression to simulations of local-type non-Gaussianity. One of 
their results is particularly interesting for our purposes: the Log-Edgeworth mass function 
is in better agreement with simulations that have a local type trispectrum with ampli- 
tude greater than (the natural size) 0(f^ L ). Since our feeder field scenario similarly has 
a boosted four-point amplitude, we compare here the expectation for the mass function 
using the Log-Edgeworth expansion. Of course, any theoretical mass function must be 
checked against and calibrated with simulations. Here we anticipate the possibility that 
the feeder type non-Gaussianity may be better predicted by the Log-Edgeworth expansion 
in analogy with results found for local-type scenarios that are designed to move away from 
the hierarchical scaling. 

The Log-Edgeworth mass function comes from truncating the expansion for Log-F, the 
logarithm of the collapse fraction, instead of F p8[|. The expression is 



n NG | _ F 

I log— Edgeworth ~ -CjXp 

no 



Fi(M) F 2 (M) _ 1 f F 1 (M) Y 

F[(M) + F>(M) F X (M)F[{M) 
F>(M) F (M)F (M) 



(5.26) 



2 



F 1 {M)+F 2 (M) / Fx(M) \ 
F (M) + \F (M)) j 

This formula works for both the hiearchical or feeder scenarios, provided the correct ex- 



pression for the Fi(M) are used; see subsection 5.3 . 

Figure |3] compares the Edgeworth and Log-Edgeworth mass functions for both scalings. 
The difference between the two expressions is much more significant for the feeder field type 
scaling. The Log-Edgeworth result at second order nearly corresponds with the fourth order 
Edgeworth result, and indicates a larger deviation between the two scalings than is evident 
at first order. 

5.5 Minkowski Functionals 

Minkowski functionals are morphological indicators of the statistical properties of a generic 
random field, /. They have been used to cross-check constraints of non-Gaussianity from 
the CMB and for weakly non-Gaussian fields can be approximated by a series in the 
moments of the distribution much like the Edgeworth expansion. General expressions for 
the non-Gaussian Minkowski functionals have been worked out in detail by Matsubara 



1 105]. CMB analysis has so far been done for Minkowski Functionals with the first order 



11 A different technique for improving the behavior of the series was proposed in |l04j ], and one might try 
to similarly extend it to non-hierarchical scalings. 
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z=0 z=0.5 



dn/dM dn/dM 




Figure 3: A comparison of the predictions from the second order Edgeworth and Log-Edgeworth 
mass functions. The ratio of the non-Gaussian to Gaussian mass functions of each type is plotted 
versus mass for redshifts z = and z = 0.5. Each panel shows the hierarchical (lower curves) 
and feeder field (upper curves) scaling assuming an equilateral type bispectrum with amplitude 
corresponding to /Jy 9 ™ = 250. The solid black lines show the second order Edgeworth result for 
hierarchical scaling. The dashed red lines are the Log-Edgeworth expansions for the hierarchical 
case and are indistinguishable from the Edgeworth for the parameters plotted here. The dashed 
black lines show the Edgeworth expansion for the feeder field scaling, and the dot-dashed blue 
lines show the Log-Edgeworth. Interestingly, the second-order Log-Edgeworth for the feeder field 
corresponds very nearly with the fourth order Edgeworth expansion (see previous figure) . 



correction [106, 107], but Matsubara has recently presented the analytic expressions needed 



for the second order terms for scenarios with hierarchical scaling |108| 1 (see also [109]). 

If / is a two-dimensional field, three Minkowski functionals for the normalized field 
v = f /a measure the fractional area of the regions above threshold (Vo(fj)), the length of 
the boundaries of regions above threshold (V\ (ft)), and the number of contours surrounding 
regions above threshold minus the number of contours surrounding regions below threshold 
(V^fj)). These quantities rely on gradients of the field, and so the expansion about the 
Gaussian will contain several statistics at each order n. For models with hierarchical 
scaling, these have been labeled as 2 variances, 3 reduced skewness parameters, 4 reduced 
kurtosis parameters, etc. : 

n = 2: a 2 = {f), o\ = -(fV 2 f) (5.27) 
„_„ q _(f) a _ (/ 2 V 2 /) q _ 2(|V/| 2 V 2 /) 



a* a 



l 

(/ 4 >c ^_(/ 3 V 2 /) c T ^ _ 2(/|V/| 2 V 2 /) c + (|V/| 4 \ 



n = 4 : K = — — , Ki = — — , K n 

Km 



It: 



a 6 ' 1 a 4 aj ' 11 a 2 af 
(|V/| 4 ) 



2a 2 af 

To reduce the notational complexity, we only give the explicit expressions relevant for the 
two-dimensional field of the CMB. In that case, there are three Minkowski functionals, 
expressed as a function of threshold Uf In a form appropriate for the hierarchical scaling, 
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the second order expressions are: 
V (v t )\hi e i 



(5.28) 
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As we did for the mass function, we can immediately rewrite the expressions for the 
Minkowski functionals in a form appropriate for non- hierarchical scaling. First, we will 
replace the reduced cumulants suitable for the hierarchical scaling by the dimensionless 
moments: 

(f) .<«=<£5™, M?> = a « v y *> (£) (5.29) 
,(/) = (/ 3 V 2 /> C U (II) = 2{/|V/| 2 V 2 /). + (IV/I 4 ). 
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Then, the second order expressions for the Minkowski functionals, appropriate for 
physics with feeder field like scaling are 



Vo(vt) | feeder — 



(5.30) 
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One might also in principal use Minkowski functionals to analyze fields other than the 
CMB [110, 111], and it is straightforward to extend our results here to those scenarios. 



6. Conclusions 

In this paper we have clarified how different classes of inflaton interactions are distinguished 
by the structure of correlation functions. We have found that couplings to other (non- 
inflaton) sectors lead to a very different structure than self- interactions. This difference 
may be relevant for LSS probes or, indeed, any observable that is sensitive to information 
about the fluctuations beyond the 3-point correlation function. Our analysis also sheds 
considerable light on the theoretical underpinnings of inflation. We have clarified the 
relation (or lack thereof) between the validity of perturbation theory and the Gaussianity 
of the PDF and uncovered novel new kinds of consistency relations among moments. Our 
results indicate that there may be interest in exploring a much wider class of non-Gaussian 
PDFs than has previously been studied. Moreover, we have computed the parametric 
scaling of n-point functions generated from inverse decay and also production of massive 
particles during inflation. 

At the technical level, we have shown how to construct the PDF for an arbitrary 
structure of cumulants and how to translate this into a prediction for the halo mass function 
and Minkowski functionals appropriate for the feeder scenario. This construction highlights 
an important point, which appears not to have been appreciated in previous literature. 
Namely, any truncated expansion of the PDF about a Gaussian ansatz makes assumptions 
about the scaling of the moments. Here we have made these assumptions explicit, showing 
how different truncations are appropriate for different microphysical inflationary scenarios. 

A detection of non-Gaussianity in the equilateral shape template would reveal inter- 
esting physics from the inflationary era. However, it could not be conclusively ascribed 
to a single effect (eg, small sound speed). We have demonstrated in this paper that a 
second possible source of nearly equilateral non-Gaussianity (which arises in a very natural 
way) differs qualitatively from small sound speed models, through the scaling of higher 
moments. 12 Differentiating between a small sound speed and feeder field scenarios would 
require at least one of three measurements: distinguishing the shapes at the level of the 
bispectrum, a measurement of the corresponding trispectrum amplitude, or a measurement 
of the tail of the halo mass function. All three are extremely challenging, but it would be 
interesting to work out which is most feasible. Finally, we point out that even if differen- 
tiating these scenarios is difficult, they demonstrate that the level of non-Gaussianity in 
the mass function or Minkowski functionals may not correlate with the amplitude of the 
three-point function in the way that has so far been assumed. Our scenarios show a wider 
range of possibilities. This analysis illustrates that great care is required to interpret what 
the statistics of rare objects imply for the physics of the primordial fluctuations. 

Our construction of a non-hierarchical PDF may have applications beyond the models 
considered in this paper. For example, one might consider a phenomenological Ansatz for 

12 Note that the detailed shape of the bispectrum and the value of the tensor-to-scalar ratio in the two 
cases may also be different JEtI]. 
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the primordial curvature perturbation of the form 



C = C 9 + <r 2 g -(a 2 g ) (6-1) 

where £ 9 and a g are statistically independent Gaussian fields. Such a phenomenological 
approach has been considered in |112| , 113 ]. (See also [114, 115] for related discussion.) 



The Ansatz fl6.1| ) is expected to give a structure of cumulants similar to the feeder mech- 
anism. However, being local in position space, Eqn. Q6.1| ) may be more easily ammenable 
to numerical simulations. 

In the body of the paper, we have emphasized that the feeder mechanism is a "single 
field" scenario because a single field sources the inflationary background and super-horizon 
fluctuations, and because there is no possibility of isocurvature. Interestingly, however, 
there are two statistically independent contributions to the fluctuations in the power spec- 



trum, as shown in Eq.( 3.24 ). This means that there may be observable stochasticity, since 
the non-Gaussianity is only correlated with a part of the two-point function (or effectively 
uncorrelated if the feeder contribution to the two-point is very small). So far, stochasticity 
has been discussed when the post-inflationary gravitational perturbation has sources from 



two different primordial fields 1 116, 117]. In the feeder field case, it is the original inflation- 
ary curvature perturbations that exhibit stochasticity, and this will carry over into LSS 
observables. Finally, we note that quasi-single field scenarios [118] may also exhibit feeder 
scaling in some regimes. 

The examples which we have studied assume that one starts with the full inflaton ac- 
tion, designed to generate both the nearly de Sitter background evolution and the spectrum 
of fluctuations. A different approach is to separate the physics of the background infla- 
tionary solution from the physics of the fluctuations - to plead ignorance about anything 
about the UV mechanism that may relate those points. This is the effective field theory of 
the fluctuations |5^], which parametrizes single-field non-Gaussianity. The usual approach 



to the effective field theory of fluctuations in single field inflation [58] (and see |119f | for 
a detailed discussion of the trispectrum) gives a hierarchical scaling of moments, by con- 
struction. There is, of course, no conflict with our results since the kinds of interactions 
that give rise to the feeder mechanism are not usually considered (although they might find 
a natural home in the framework established in |120| ]). One might imagine trying to incor- 
porate analogous effects into the "single field" framework by replacing the feeder field with 
an external source. Such a term would give rise to a tadpole contribution that modifies the 
Friedman equation and is usually omitted. In the feeder scenarios that we have considered, 
one can work in a consistent regime where modifications to the Friedman equation are 
negligible, although there is still a non-trivial impact on cosmological fluctuations. 

In this work we have considered only two realizations of the feeder mechanism, arguing 
that both give rise to a similar structure of correlation functions. We expect that this 
structure is quite general. It would be interesting to construct more explicit models of this 
type and verify this expectation. 
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A. Taxonomy of Inflation Models 

A.l A First Taste of Inflation: Vanilla Flavoured Models 

In the vanilla scenario, the inflaton is effectively decoupled from all other species of particles. 
Self interactions of the scalar curvature mode arise either through the inflationary potential 
or through gravitational couplings. In both cases, the interaction strength is controlled by 
slow roll parameters, leading to an important consistency relation between the size of non- 
Gaussianity and the background dynamics. That is, the cosmological fluctuations in the 
vanilla scenario deviate from Gaussianity to roughly the same extent that the background 
deviates from pure de Sitter §(| 0> Hi- 



The vanilla model (2.1) is simple, but not entirely satisfying for several reasons. The 
first reason is phenomenological. The observational signatures of the vanilla model are 
minimal: we have already observed the the amplitude of the curvature fluctuations (the 
ratio of the scales H and yfeM p ) and have perhaps obtained some evidence for a red spectral 
index. The only observable left is the amplitude of gravitational waves (which determines 
H); however, this may well be too small to measure in the next decade. Primordial non- 
Gaussianity, although present, will be swamped by post-inflationary gravitational evolution 



1 121, 122 , [L23|, 124, 125] and so is also extremely difficult to measure in the immediate 



future. Although the vanilla model is consistent with the data IS], it will be very hard to 
say how it can be genuinely tested in the scalar sector. On the other hand, extensions of 
the vanilla scenario that give observably large correlations beyond the two-point are much 
more falsifiable. 

The second reason for exploring alternatives to the vanilla scenario is theoretical. If 
the basic framework of the inflationary paradigm is correct, it may perhaps be surprising to 
find only one degree of freedom, decoupled from the rest of particle physics, and no relevant 
energy scales in the problem besides those associated with gravity or the slow roll potential. 
Indeed, efforts to embed inflation within a realistic particle physics framework (such as 
string theory, the MSSM or supergravity) typically contain many dynamical scalars, in 
addition to various couplings between the inflationary sector and matter fields. (We use 
the word "matter" to denote any field which is not involved in driving the inflationary 
expansion.) In section |2] we provided a simple and explicit example illustrating how a 
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rich, non-vanilla interaction structure arises automatically once we try to realize single- 
field inflation in a more complete particle physics framework. These "new" interactions 
are controlled by a scale, / <C M p , and they can lead to observationally interesting non- 
Gaussianity. 

Once we move beyond the vanilla scenario, there are clearly several possibilities: the 
background evolution and/or the curvature fluctuations on superhorizon scales may have 
contributions from multiple fields, and the fields sourcing the background may or may not 
be the same as those sourcing the curvature fluctuations. 

A. 2 Classifying Models: I am 32 flavours and then some 

We are now in a position to provide a taxonomy of known inflationary models. To avoid 
confusion, we remind the reader that our definition of "single field" still allows a "single 
field" inflaton to couple directly to matter fields. 



1. Decoupled Single Field Inflation: There is only a single scalar degree of freedom 
which sources both the background evolution and also the curvature fluctuations. The 
inflaton is effectively decoupled from matter fields during inflation. Observable non- 
Gaussianity may be generated if there are self-interactions - beyond the usual slow- 
roll/gravitationally suppressed "vanilla" couplings - controlled by a scale / <C M p . 
In this case, the bispectrum is closer to equilateral than local, since it is generated by 
physical processes on scales of order the horizon or smaller. There are no iso-curvature 
perturbations and no super- horizon evolution of £. 

2. Coupled Single Field Inflation: There is a single scalar inflaton degree of freedom 
and the background dynamics are of the slow roll type. Nontrivial couplings to 
matter fields may generate significant non-Gaussianity via the feeder mechanism; see 



subsection 2.5. The bispectrum is closer to the equilateral than local shape, since 
the relevant interactions are most significant near horizon crossing. There are no 
appreciable large-scale iso-curvature perturbations or super-horizon evolution. 

Dissipative Models: There is a single scalar inflaton degree of freedom which is 
very strongly coupled to matter fields. The production of matter quanta is so strong 
that dissipative effects contribute the dominant source of friction for the motion 



of the inflaton. Examples include trapped inflation [126, [n] (see also [gg, 87] for 
generalizations) and the Anber Sz Sorbo model [|35j| , both of which may be thought 
of as an extreme limit of the feeder mechanism. Non-Gaussianity is expected to be 
closer to equilateral than local. (Warm inflation [ |127| ] might also be included in this 
category.) 

4. Curvaton: Here the link between the background evolution and the curvature fluc- 
tuations is loosened. The fluctuations of a second field, irrelevant during inflation, 
generate the dominant contribution to the curvature perturbations at the end of in- 
flation. Non-Gaussianity is of the local type and isocurvature modes may be relevant 
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depending on details of the scenario. When the curvaton does not completely dom- 
inate the fluctuation spectrum (a mixed inflaton/curvaton model) the perturbation 
spectrum can be fundamentally multi-field. 

5. Multi-field inflation: Here more than one field contributes both to the background 
evolution and to the fluctuations. Non-Gaussianity is generated by super-horizon 
evolution and so is of the local type. Isocurvature modes may be relevant depending 
on details of the scenario. 

See also Table The categories which we have enumerated here are not meant to 
exhaust all possibilities. Rather, these fundamental mechanisms should be understood as 
basic building blocks which can be assembled in various combinations to construct ever 
more elaborate models. 

Table 1: Classifying inflationary scenarios 



Model 


Background 


Super-horizon ( 


Non- Gauss. 


Iso-Curv. 


Single field? 


vanilla 


single field 


same field 


not local 


none 


V 


feeder 


single field 


same field 


not local 


negligible 


V 


dissipative 


single field 


same field 


not local (?) 


negligible (?) 


V 


curvaton 


single field 


different field 


local 


yes 


X 


multi-field 


many fields 


many fields 


local 


yes 


X 



B. Non-Gaussian Correlators from Inverse Decay 



In subsection 2T we introduced a general effective field theory description (23) which we 



expect to be valid in any model of single-field inflation that is characterized by an underlying 



shift symmetry. In subsection we discussed three distinct mechanisms which operate 
in this theory, each of which may give rise to observationally interesting non-Gaussianity. 
The first two mechanisms - sound speed and resonance effects - have already received 
considerable attention in the literature. The third mechanism, on the other hand, is newer 



1 44]. In this appendix we compute the parametric scaling of all non-Gaussian n-point 



functions generated by this mechanism. The calculation is a straightforward application of 



the formalism detailed in [57| and the reader is referred to that paper for more details. 



B.l Cosmo logical Perturbation Theory 

Following p4|], we focus on the theory ( |2.9| ) in the case where both higher derivative cor- 
rections and also resonance effects are negligible. We adopt the flat slicing and Coulomb 
gauge. By integrating out the shift and lapse functions one may derive an action for the 
perturbations 5ip(t,x.) = (p(t,x) — </>(£) and A^(t,x). Keeping terms up to third order, this 
may be written as |57[| : 

S = S + SS grav (B-l) 
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The leading order action contains the kinetic terms for the gauge field and inflaton, in 
addition to the pseudo-scalar coupling: 



So 



+ 



drdrx 



drd x 



<5<S 



'\2 



a 2 M< 



SpeijkA'idjAk 



t.A\A\ — -diAjdiAj 



(B-2) 



° Um p 



and d 2 is the inverse Laplacian satisfying 



where we have defined M 2 = V" 
dr*di&f = f. 

The full third-order action ( B-l| ) also contains cubic interactions which are associated 
with the nonlinearity of gravity and hence are suppressed, both by the Planck scale and 
also by slow roll parameters. In (B-l) we have denoted these by 5S grijiV . Explicitly, these 



are |28], y|: 



SS, 



grav 



+ 



drd 6 x 



drd 6 x 



M 



5.p 



v 
2e 



-\A\A\ - X -F l3 F l3 + X -d- 2 d T d^F l3 A))^ 



-5ip'did- 2 (5ip')di5<p - —5<p(5ip' 



f\2 



-S^(V5^) 2 (B-3) 



Evidently, the strength of the interactions encoded in ( |B-3| ) is set by the coupling constant 
This should be compared to the scale j which controls the strength of the pseudo- 



scalar interaction in flB-2|) . As discussed in section g, the natural expectation from effective 
field theory is a = 0(1) and / -C M p . In this case, we have j ^> and the gravitational 

couplings 55 g rav may be neglected. Quantitatively, we expect e 1//2 ~ 10 _1 for a large- 
field inflation model, whereas aM p /f ~ 10 2 whenever non-Gaussianity is observationally 
interesting [H4J|. Hence, we expect that the gravitational couplings are negligible in the 
most interesting scenario. 

Notice that the terms on the second line of (|B-3| ) are just the usual trilinear scalar 
interactions which would arising in ordinary slow roll inflation, even in the absence of the 
gauge field. These are known to contribute Af NL = C(e, 77) ||, 13' H> B HI- It follows 
that such interactions may be neglected whenever /jvx ~ 0(1)- 

B.2 Production of Gauge Fluctuations 

The key aspect of the dynamics of the model ( |B-1| ) is the production of fluctuations, SA^, 
due to the homogeneous motion of the inflaton. To capture this effect we study the equation 
of motion for the gauge field fluctuations in the homogeneous background 4>(t): 



A" 



V 2 A 



f 



■Vxvl = 



(B-4) 



The last term arises due to the coupling of the gauge fluctuations to the time-dependent 
condensate. This term breaks the conformal invariance of the gauge field sector and, 
roughly speaking, becomes important for wave- numbers k <, a<j)' //. 
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We decompose the q-number field A(t, x) as 



^3 [e A (k)a A (k)A A (r, k)e ik ' x + h.c.] (B-5) 

where "h.c." denotes the Hermitian conjugate of the preceding term, the annihilation/creation 
operators obey 

a A (k),aUk')l = (2^) 3 W (3) (k - k') (B-6) 



and e A are circular polarization vectors (see |57|] for more details). Using the decomposition 
QB-5| ) the equation of motion for the c-number mode functions in the classical background 
(p(t) takes the form 



9 - + k* ± ^ 



A ± (r,k) = 0, t = ^fe (B-7) 



During inflation the parameter £ may be treated as constant, as its time variation is 
subleading in a slow roll expansion. 

From equation ( |B-7| ) we see that one of the polarization states experiences a tachyonic 
instability for k/(aH) < 2£. (Notice that this effect could arise also in flat space, the depen- 
dence of £ on the expansion rate H arises due through time dependence of 4>(t).) Without 
loss of generality, we assume that <j) > during inflation, so that the mode exhibiting the 
instability is A + . The correctly normalized solutions of (|B-7|) are well approximated by 



in the interval (8£) _1 < k/(aH) < 2£ of phase space that accounts for most of the power 
in the produced gauge fluctuations. The phase space of growing modes is non-vanishing for 
£ }t 0(1), which we assume throughout. Notice the exponential enhancement e n ^ in the 
solution ( |B-§| ). This arises due to the tachyonic instability and reflects significant nonper- 
turbative gauge particle production in the regime £ > 1. On the other hand, the production 
of gauge field fluctuations is uninterestingly small for £ < 1. The other polarization state, 
A-(r, k), is not produced and can therefore be ignored. 

B.3 Inverse Decay Effects 

In the last subsection, we showed that the homogeneous motion of the inflaton leads to 
a tachyonic production of gauge fluctuations. Throughout this paper, we work in the 
conventional slow roll regime where the backreaction of these produced fluctuations on the 
homogeneous background is negligible; more on this later. 13 Nevertheless, the production 
of gauge quanta still has an important effect on the observable cosmological fluctuations in 
the model: the produced gauge fluctuations may source inflaton perturbations via inverse 
decay, 5 A + 5 A — > 5ip. 



3 See |B5[| for an alternative scenario where backreaction effects are important. 
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We decompose the q-number inflaton fluctuation as 



5<p(r, x) 



d 3 k Q k (r) 



„ik-x 



(2tt) 3 o(t) 

The equation of motion for the Fourier modes is 



^ + k 2 + a 2 M 2 



a" 



J k (r) = a s 



4/ 



a 

j3„„— tkx 



drxe~ 



Qk(r) = J k (r) 



(B-9) 

(B-10) 
(B-ll) 



where the quantity Fi 7 is constructed from the solutions ( |B-§| ) . The solution of ( B-1C] ) 
may be split into two parts 



homogeneous particular 



(B-12) 



The homogeneous solution corresponds, physically, to the usual vacuum fluctuations from 
inflation. The particular solution, on the other hand, may be interpreted as arising due to 
inverse decay processes. 

The usual vacuum fluctuations from inflation are Gaussian, to a very good approxi- 
mation 2C, 21, [22]. Hence, the homogeneous term may be expanded as 

gjT(r) = b(k)u k (r) + 6t(-kK(r) (B-13) 

where the c-number modes u k are the well-known homogeneous solutions of flBTLOP . The 
inflaton ladder operators obey 



6(k),6t(k / ) = (2vr) 3 ^ 3 )(k-k') 
and commute with the ladder operators of the gauge field 



[6(k),a A (k')] = &(k),4( k " 



(B-14) 



(B-15) 



Next, we turn our attention to the particular solution of ( |B-10| ), which may be written 



as 



/^inv.dec f_\ 



f dr'G k (T,T')MT') 

J —00 



(B-16) 



G k (T,T') = 5(r-r') 



(B-17) 



where the source term was defined in ( |B-11 ) and G k (r, r') is the retarded Green function 
which solves 

% + + — 

a 

The n-point correlation functions take the form 

(Q inv.dec (T) ... Q inv.dec (r)) = 

X (J kl (Tl)---Jk n (Tn)) (B-18) 



Y[ / dTiG ki (T,Ti 
.1=1 
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Ultimately, we are interesting in the co-moving curvature perturbation, £(r, x), which 
is related to the field fluctuation as ( |2.4D . Working in Fourier space, the contribution to £k 
from inverse decay effects is 

C mv.dec (r) = _#Qinv.dec( T ) (B _ lg) 



Following closely the formalism detailed in [57], it can be shown that, on large scales 



-kr <C 1, this may be written in a suggestive form: 

cr- ( .) - 0.0031 ^euf^^ 



2vr|0| / V e 4 J ^ J (2vr)3 
a+(q)a+(k - q) + a f + (-q)a + (k - q) + a\. (-k + q)a+(q) (B-20) 

+a f + (-q)aU-k + q) 

Here we have introduced dimensionless ratios k = k/k, q* = q/k and also defined the 
quantity 

\n |V4|& - n I 1 / 4 

&) = 4.2 • 10 2 ei (&) ei (fc - &) , |g *' 1 ^ -j (B-21) 

|£|l/2 + |£_£|l/2 



which is dimensionless and has no dependence of model parameters whatsoever. We have 
normalized the ^-function so that 



/ 



< B - 22 » 



which cancels the combinatorial factor arising from the Wick contractions in the 2-point 
correlator. 

Equation ( B-20) ) has been written in a way which makes manifest the dependence of 
££ lv " dec (r) on model parameters. The prefactor multiplying ( p-20 ) can be understood as 



follows. The factor which would give the amplitude of fluctuations in the vanilla 

2n\(p\ 

scenario, can be traced to the relation (|2.4|). The quantity ^j-, on the other hand, is the 
natural dimensionless measure of the strength of the pseudo-scalar interaction that gives 
rise to inverse decay effects. Finally, the contribution is associated with the time 
dependence of the gauge modes ( |B-8| ). 14 This last term depends on <fi, H and a//, but 
only through the dimensionless combination £ that controls the effective mass of the gauge 
field; see equation flB-7| ). 

Two things are evident from equation ( |B-20 ). First, the particular solution ( B-16| ) is 



statistically independent of the homogeneous solution ( B-13| ). Indeed, the former can be 
expanded in terms of the annihilation/creation operators a^(k),at(k) associated with the 
gauge field while the latter is expanded in terms of the annihilation/creation operators 



In order to derive (B-20) we work in the limit £ 2> 1. At smaller values of £ the power-law dependence 
of this term will be modified, as will the factors order unity. Neither of these modifications will have any 
impact on our qualitative results for the structure of the correlation functions. 
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6(k), &t(k) associated with the inflaton vacuum fluctuations. As we pointed out, these two 
sets of operators commute with one another. 

The second notable feature of equation ( |B-20| ) is its non-Gaussianity, which is evident 
from the quadratic dependence on the annihilation/creation operators a,\(k), at (k). This 
is easily understood: the particular solution Q™- dcc is bi-linear in the Gaussian field SA^ 
and hence we expect this to behave, heuristically, as the square of a Gaussian field. 

B.4 Parametric Scaling of the Correlation Functions 

We now turn our attention to the n-point functions of the co-moving curvature perturbation 



H 



Ck(r) = -— Qr(r) + QZ 
ac 



inv.dec . 



r) 



(B-23) 



Before considering non-Gaussian correlators, the 2-point function deserves special atten- 



tion. We have (44 



(CkCk 
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n 2 /A2 



H 4 



//^>vac /^ivac\ , //^inv.dec /-Mnv.dcc\ 

Wk Vk' / + Wk Vk' I 



(2tt)2|^|5 



2vr 2 + 



0.0031 • 10" 



{ aH e 2 < 



k 3 



6® (k + k') (B-24) 



Notice that there is no "cross-term" because the two contributions to £k are statistically 
independent. The first term in the square braces comes from the quantum vacuum fluctu- 
ations, while the second term arises due to inverse decay effects. 

We now turn our attention to a parametric computation of all connected n-point 
functions, with n > 3. These correlators encode departures from Gaussianity, whereas 
the usual vacuum fluctuations from slow roll inflation are known to be Gaussian to high 
accuracy [^, |2l], 22] . Hence, it will be a good approximation to consider only the particular 
solution of ( |B-10D and neglect the homogeneous solution. In this approximation, and in 
the limit -fcr < 1. we can write 



(CkiCka • • "CO 



0.0031- 



H 2 aH e 2 ^ 



2n\ 



f e 4 
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X (aqiOki- 



qi 



-k n +q n 



)B-25) 



where we have employed equation ( |B-20| ) . The expectation value on the last line involves 
a sum of products of In annihilation/creation operators. On evaluating the Wick contrac- 
tions, these give n delta functions, each with a factor of (2-7r) 3 . By translational invariance, 
one of these must be (2vr) 3 ^ 3 ) ki). The remaining n — 1 delta functions kill off all but 
one of the integrals J . The final integral over internal momenta gives a factor when 



put into dimensionless form: f J^L 



k 3 J 



Thus, the final correlator must scale like 



<(an 



0.0031- 



H 2 aHe 2 ^ 



2vr| 



/ e 4 



(2*)W J> 



(B-26) 
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The factor fc -3 ^" 1 ) simply reflects the scale invariance of the produced fluctuations. The 
coefficients A n are dimensionless numbers, independent of model parameters, which encode 
the integration over ^-functions and also combinatorial factors from the Wick contractions. 
Our normalization of the -(/'-function has been chosen so that A2 = 1 and it may be checked 
explicitly that A3 ~ 2.2. We expect that A n = 0(1) for higher order n also. Note that 
at very large values of n, combinatorial factors may become important and influence the 
scaling. This issue, which is by no means special to the inverse decay scenario, is left to 
future studies. 

For later convenience, we write our key results in configuration space: The variance is 
given by 



(C 2 (x)> 



The cumulants are given by 



dk 



H 4 



(2tt)2|# 



1 + 



1 
2^2 



0.0031- 
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(B-27) 
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(2vr) 3 ^ n 0.0031 
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(B-28) 



B.5 Computational Control 



A number of simplifying assumptions were made in our derivation of the key result ( B-26| ) 
for the parametric scaling of the non-Gaussian correlation functions from inverse decay 
effects. In this subsection, we justify those approximations in more detail and draw atten- 
tion to the various constraints on model parameters that are imposed by the requirement 
of computational control. 

As discussed in the last section, we require H <C / <C M p in order for the effective 
field theory description ( |B-1| ) to be valid. In addition to this, we have made two main 
simplifying assumptions in deriving ( B-26| ). Namely: 



1. Neglect of Gravitational Interactions: We have neglected the gravitational in- 
teractions encoded in 5S gTStv ; see equation ( p3-3| ) . 

2. Neglect of Backreaction: We have neglected the backreaction of the produced 
gauge fluctuations on the homogeneous dynamics of 4>(t) and a(t). 

1/2 

We have already argued that the former assumption is justified whenever j 3> j^-; see the 
discussion below flB-3j ). The latter condition, on the other hand, requires somewhat more 
care. 

There are two distinct backreaction effects. First, the gauge fluctuations SA^ are 
produced at the expense of the kinetic energy of the inflaton condensate, 4>{t). This effect 
is encoded in the mean field equation 



+ 3H4> + V'(4>) = j{E-B) 



(B-29) 
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where E = — -^jA' and B = 4yV x A are the physical electric and magnetic components of 
A^. In order to trust the standard slow roll result 3H<j) = — V we require \ V'\ ^> j(E ■ B). 
This condition may be shown to be equivalent to [57] 

H 2 aHe 2 ^ 



0.0031 



2vr|0| / £ 4 



< 0(1) 



(B-30) 



The second backreaction effect arises because the produced gauge fluctuations con- 
tribute to the total energy density of the universe, modifying the Friedman equation: 



3H 2 



1 

Ml 



^ 2 + V( ( p) + 1 -(E 2 +B 2 ) 



(B-31) 



It can be shown that the energy density in gauge field fluctuations is negligible, (E 2 +B 2 ) -C 
V, whenever the condition ( |B-30 ) is satisfied. 

Notice that the backreaction bound ( B-30 ) constraints the same parameter which 
controls the size of the non-Gaussian n-point correlation functions; see equation ( B-26| ). 
Indeed, this condition could be re-phrased as 



inv.dec 



,V2 « ! 



(B-32) 



This is not surprising since, after all, we expect backreaction effects to be negligible when 
fluctuations are small. 

Our solution ( |B-g| ) for the produced gauge fluctuations was obtained neglecting in- 
homogeneities of ip. As a final consistency check, let us verify that this approximation 
is valid, even in the inverse decay dominated regime. The equations of motion of the 
electromagnetic field can be written in the form 



E' + 2aHE - V x B = -jv'B ~ J^p x E 



V • E 



a /, 



An interative solution of the equations of motion will converge provided 

1/2 



^( J B 2 ) 1/2 »y((^x/: 




(B-33) 
(B-34) 

(B-35) 



We can estimate 
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e(c 2 ) 1/2 #<£ 2 > 1/2 



(B-36) 



Hence, we arrive at the condition (C 2 ) 1 ^ 2 <C which is easily satisfied whenever back- 
reaction effects are small (recall that £ •> O(l) in the most interesting region of paramter 
space). 
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C. Non-Gaussian Correlators from Rescattering 



In this appendix we estimate the cumulants for the model (4.1). Decomposing the q- number 
inflaton fluctuations as QB-9| ) the equation of motion (4.3) takes the form 



dl + k 2 + 



a m 



a 
a 



Qk(r) = J fc (r) (C-l) 
J k (r) = -a 3 gk 2 Ar) f d 3 xe~^ [ X \r, x) - ( X \r, x))] (C-2) 

where k+ = \J g\(j)\, t(r) = j T a(r')dr / and we have arbitrarily set the origin of time so that 
4>{t) = (j>Q at t = 0. The solution of ( C-l ) may be written as 

Qk(r)= QT(t) +QTM (C-3) 

homogeneous particular 

The homogeneous solution corresponds, physically, to the usual vacuum fluctuations from 
inflation. This is given explicitly by ( |B-13| ). The particular solution, on the other hand, 
may be interpreted as arising due to rescattering. 

The vacuum fluctuations are Gaussian to very high accuracy, hence, to compute the 
cumulants with n > 3 we can focus only on the particular solution of ( p-l| ) . The correlation 
function of the curvature perturbation is 



(Ck!Ck 2 •••Ck n (r)) c 
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x (J kl {n)Jk 2 (T2) ■ ■ ■ JkM) (C-4) 



Ref. 1 65] provides a detailed account of x-P ar ticl e production in an expanding universe, 
showing explicitly how to solve for the produced fluctuations, x( T > x )- For our purposes, it 
is sufficient to note that the n-point correlator of the source term can be estimated as 

(JkAh^kM ■ ■ ■ Jk n (tn)} ~ g n kle-* ck2 l^ x JJefe) x 8® {^2 k ^j ( C " 5 ) 

where c is some order unity number. The product over step functions enforces the fact 
that the source Jki(U) turns on only for t{ > 0, after particle production has occured. 
The estimate ( p-5[ ) would recieve order unity corrections in a more careful computation, 
however, this expression correctly captures the parametric dependence and will be quite 
sufficient for our purposes. 

Using the estimate (|C-5| ) the n-point correlation function becomes 
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(C-6) 
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where tq is the value of the conformal time variable when (p = <p$. If we set a = 1 at t = 
then To = —1/H. The time dependence arises through integrals (called I2 in |35|]) of the 
form 

/ dn 

J-l/H 



G k (T,n) J_ f (Jl 



a(r) 



-f , 

H 2 J \H 



(C-7) 



where we take the limit — kr — > and introduce the function 

fix 



X 



1 f x dx' 1 
-g j — - [sin(x') — x'cos(x / )] = — [Si(x) — sin(x) 







(C-8) 



Here Si(x) is the sine integral. 

Putting everything together, we estimate the momentum space correlator as 
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As expected, the fluctuations from inverse decay are not scale invariant, rather, the effect 
peaks near scales leaving the horizon when particle production occured. 
The cumulants are 
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Most of the support for this integral comes from scales k ~ clqH, corresponding to the 
location of the non-Gaussian feature. (Recall that we have set ao = a(t = 0) = 1.) We can 
smooth the fluctuations on the scale of the feature simply by constraining the domain of 
integration to lie close to the peak of the integrand. Proceeding in this way gives 



(C(x)) Rb ~ a n ^ 



(2vr<7) 



H 2 



2vr|0| 



(C-ll) 



where the coefficients a n are independent of model parameters and Rb ~ H . The variance 
recieves important contributions both from rescattering effects and also from the usual 
vacuum fluctuations. Summing up both contributions, we have 

k 3 



<a*)>* 



(C-12) 



(2^) 2 2 

Both equations (|C-11|) and (|C-12| ) are very similar to the analogous result for inverse decay 
effects. This is not surprising, given the qualitative similarity of the underlying physics. 
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